Spark Ignition of Monodisperse Fuel Sprays 


by 


Allen M. Danis 


Submitted in partial fulfillment of the requirements for 

the degree of 


Doctor of Philosophy 


Approved by 


■n j-x, £ 




Dr. Nicholas P. Cernansky a/rd Dr. Izak Namer 
Co-Principal Investigators 


Research Sponsored by the NASA Lewis Research Center 
Grant No. NAG 3-382 and by Drexel University 


Department of Mechanical Engineering and Mechanics 
Drexel University, Philadelphia, PA 19104 


September 1987 



ABSTRACT 


A study of spark ignition energy requirements was conducted with a 
monodisperse spray system allowing independent control of droplet size, 
equivalence ratio and fuel type. Minimum ignition energies were measured 
for n-heptane and methanol sprays characterized at the spark gap in terms of 
droplet diameter, equivalence ratio (number density) and extent of 
prevaporization. The droplet diameter was varied between 30 and 57 ym, 
and the equivalence ratios ranged from 0.45 to 1.8. The extent of 

prevaporization of the sprays was held approximately constant for each fuel, 
0.45 for n-heptane and 0.34 for methanol. Droplet Reynolds numbers of the 
sprays ranged from 0.1 to 1.0, and interdroplet spacings varied between 15 
and 30 droplet diameters. In addition to sprays, minimum ignition energies 
were measured for completely prevaporized mixtures of the same fuels over 
a range of equivalence ratios to provide data at the lower limit of droplet 
size. 

Results showed that spray ignition was enhanced with decreasing droplet 
size and increasing equivalence ratio over the ranges of these parameters 
studied. No optimum equivalence ratio for ignition was observed for sprays, 
but an optimum did occur at 0 * 1 .8 for prevaporfzed n-heptane ignition. A 
corresponding optimum for prevaporized methanol could not be obtained due 
to excessive fuel condensation. 

By comparing spray and prevaporized ignition results, the existence of an 
optimum droplet size for ignition (below 30 pm) was indicated for both fuels. 
This optimum was attributed to the maximum flame propagation rates 


observed by previous researchers for sprays In this droplet size range ( 10 - 
40 wm), where both premixed and diffusive types of combustion occur 
simultaneously. Extension of the lean Ignition limits was also observed for 
sprays compared to prevaporized mixtures. This extension occured for both 
fuels over the range of droplet sizes studied, and was attributed to fuel 
distribution effects (i.e. fuel rich regions surrounding evaporating 
droplets). Both the optimum droplet size and the extension of the lean 
ignition limit indicate the importance of the gas-phase stoichiometry in the 
interdroplet spacing to spray ignition. 

Fuel volatility was seen to be a critical factor in spray ignition. The 
more volatile n-heptane sprays required roughly one-fifth the ignition energy 
required by methanol sprays of corresponding size and equivalence ratio. 
For prevaporized ignition, however, activation energy was the important fuel 
property. Prevaporized methanol, with its lower activation energy, required 
roughly 50* less ignition energy than prevaporized n-heptane. 

The spray ignition results were analyzed using two different empirical 
ignition models for quiescent mixtures. Both models employed a 
characteristic time approach, relating the time required for fuel evaporation 
in the spark kernel to the quenching time of the kernel by thermal conduction 
to its surroundings. The major difference between the models was that while 
the general ignition model developed by Ballal and Lefebvre equated the two 
times, the Characteristic Time Model of Peters and Mellor stated that they 
were only proportional and determined the correlation between them. Both 
models accurately predicted the experimental ignition energies for the 
majority of the spray conditions. Their performance deteriorated for lean 
equivalence ratios (< 0.7) and smaller droplet sizes (< 40 ym) however. 
This was probably because the models did not account for the time required 



for chemical reaction In the spark kernel, which becomes Important for 
leaner and smaller droplet size sprays. 

Spray ignition was observed to be probabilistic in nature, and ignition was 
quantified in terms of an ignition frequency for a given spark energy. The 
criterion for the minimum Ignition energy was the spark energy which 
produced an ignition frequency of 50*. A model was developed to predict 
ignition frequencies based on the variation in spark energy and equivalence 
ratio (number density) in the spark gap. Random normal distributions of 
these two parameters were generated with a Monte-Carlo routine and used in 
conjunction with the Characteristic Time Model for ignition. The resulting 
ignition frequency simulations were nearly identical to the experimentally 
observed values, indicating that the probabilistic nature of spray ignition is 
the result of variations in the energy levels of individual sparks and the 
random distribution of droplets in a fuel spray. 
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ABSTRACT 


SPARK IGNITION OF MONODISPERSE FUEL SPRAYS 

Alien M. Danis 
Nicholas P. Cernansky 
Izak Namer 

A study of spark ignition energy requirements was conducted with a 
monodisperse spray system allowing independent control of droplet size, 
equivalence ratio and fuel type. Minimum ignition energies were measured 
for n-heptane and methanol sprays characterized at the spark gap in terms of 
droplet diameter, equivalence ratio (number density) and extent of 
prevaporization. The droplet diameter was varied between 30 and 57 pm, 
and the equivalence ratios ranged from 0.45 to 1.8. The extent of 
prevaporization of the sprays was held approximately constant for each fuel, 
0.45 for n-heptane and 0.34 for methanol. Droplet Reynolds numbers of the 
sprays ranged from 0.1 to 1.0, and interdroplet spacings varied between 15 
and 30 droplet diameters. In addition to sprays, minimum ignition energies 
were measured for completely prevaporized mixtures of the same fuels over 
a range of equivalence ratios to provide data at the lower limit of droplet 
size. 

Results showed that spray ignition was enhanced with decreasing droplet 
size and increasing equivalence ratio over the ranges of these parameters 
studied. No optimum equivalence ratio for ignition was observed for sprays, 



but an optimum did occur at 0 * 1.8 for prevaporized n-heptane Ignition. A 
corresponding optimum for prevaporized methanol could not be obtained due 
to excessive fuel condensation. 

By comparing spray and prevaporized ignition results, the existence of an 
optimum droplet size for ignition (below 30 pm) was indicated for both fuels. 
This optimum was attributed to the maximum flame propagation rates 
observed by previous researchers for sprays in this droplet size range ( 10 - 
40 wm), where both premixed and diffusive types of combustion occur 
simultaneously. Extension of the lean ignition limits was also observed for 
sprays compared to prevaporized mixtures. This extension occured for both 
fuels over the range of droplet sizes studied, and was attributed to fuel 
distribution effects (i.e. fuel rich regions surrounding evaporating 
droplets). Both the optimum droplet size and the extension of the lean 
ignition limit indicate the importance of the gas-phase stoichiometry in thq 
interdroplet spacing to spray ignition. 

Fuel volatility was seen to be a critical factor in spray Ignition. The 
more volatile n-heptane sprays required roughly one-fifth the ignition energy 
required by methanol sprays of corresponding size and equivalence ratio. 
For prevaporized ignition, however, activation energy was the important fuel 
property. Prevaporized methanol, with its lower activation energy, required 
roughly 50% less ignition energy than prevaporized n-heptane. 

The spray ignition results were analyzed using two different empirical 
ignition models for quiescent mixtures. Both models employed a 
characteristic time approach, relating the time required for fuel evaporation 
in the spark kernel to the quenching time of the kernel by thermal conduction 
to its surroundings. The major difference between the models was that while 
the general ignition model developed by Ballal and Lefebvre equated the two 
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times, the Characteristic Time Model of Peters and Mellor stated that they 
were only proportional and determined the correlation between them. Both 
models accurately predicted the experimental ignition energies for the 
majority of the spray conditions. Their performance deteriorated for lean 
equivalence ratios (< 0.7) and smaller droplet sizes (< 40 Jim) however. 
This was probably because the models did not account for the time required 
for chemical reaction in the spark kernel, which becomes important for 
leaner and smaller droplet size sprays. 

Spray ignition was observed to be probabilistic in nature, and ignition was 
quantified in terms of an ignition frequency for a given spark energy. The 
criterion for the minimum ignition energy was the spark energy which 
produced an ignition frequency of 50*. A model was developed to predict 
ignition frequencies based on the variation in spark energy and equivalence 
ratio (number density) in the spark gap. Random normal distributions of 
these two parameters were generated with a Monte-Carlo routine and used in 
conjunction with the Characteristic Time Model for ignition. The resulting 
ignition frequency simulations were nearly identical to the experimentally 
observed values, indicating that the probabilistic nature of spray ignition is 
the result of variations in the energy levels of individual sparks and the 
random distribution of droplets in a fuel spray. 



CHAPTER 1 


INTRODUCTION 

The combustion of liquid hydrocarbon fuels is the principal source of 
energy production in today’s society. Although their percentage of total 
worldwide energy consumption dropped from 48% in 1975 to 40% in 1985, liquid 
petroleum products still account for the major portion of energy consumed in 
the U.S. and worldwide (DOE/EIA, 1985). Roughly 50% of these liquid 
hydrocarbon fuels are burned in spray form, by combustion devices such as 
furnaces and gas turbine and diesel engines. 

Studying the ignition characteristics of liquid hydrocarbon sprays is 
important for several reasons. Since many of the combustion devices 
mentioned above use spark ignition systems, the amount of spark energy 
required for successful ignition over a range of operating conditions must be 
known. Information about the ignition characteristics of fuel sprays and 
vapor is also useful for the development of fundamental models of ignition, 
combustion and flame propagation. Finally, knowledge of ignition 
mechanisms is necessary for the prevention of unwanted explosions and 
fires. 

A local source of ignition, such as an electric spark, initiates chemical 
reaction by energizing a small volume (or kernel) of combustible mixture 
around the spark. Energy is lost from this spark kernel by heat transfer and 
transport processes, and the chemical energy of the mixture in the kernel 
must be converted into thermal energy at a sufficient rate to overcome these 


losses. If this occurs, the kernel will grow into a steadily propagating flame 
and ignition is successful. 

The many interacting factors which influence spark ignition of fuel sprays 
can be grouped into the three basic categories consisting of spray 
parameters, physical parameters and spark parameters. Spray parameters 
include the size or size distribution of droplets in the spray, the interdroplet 
spacing, the amount of evaporated fuel and oxidizer present in the 
interdroplet spacing, and the fuel properties. Physical parameters which 
can affect ignition are flow characteristics such as the spray velocity and 
level of turbulence, as well as the local temperature, density and pressure 
at the point of ignition. Spark parameters such as the energy and duration of 
the spark and the electrode spacing and configuration relative to the spray 
must also be considered. 

Extensive studies have been undertaken previously to measure the 
minimum spark ignition energy required for a variety of homogeneous 
fuel-air mixtures. However, measuring the ignition energy of heterogeneous 
fuel-air mixtures is a much more complex task due to the experimental 
difficulties in controlling and quantifying the spray parameters mentioned 
above. The ignition of polydisperse, or multiple droplet size sprays has been 
studied for a number of fuels, but only over limited ranges of equivalence 
ratio, ignition data for monodisperse, or single droplet size sprays is even 
more sparse, with studies covering only one fuel and a limited size range. 
This study was undertaken to provide ignition energy data for monodisperse 
sprays over a range of drop sizes and equivalence ratios more typical of 
those found in practical combustors. 

A well characterized spray facility was used to simulate the spark 
ignition process which occurs in practical combustors, such as gas turbine 


engines. The major experimental variables included fuel type, droplet size 
and stoichiometry. The main experimental measurements made in this study 
were ignition frequency and spark energy (from which ignition energy was 
determined), droplet number density (from which the stoichiometry was 
determined) and droplet size. Ignition energy measurements were also 
performed for prevaporized, premtxed cases to provide the lower bound of 
zero droplet size. The analytical portion of the study involved predicting the 
minimum ignition energy of the fuel sprays as a function of fuel type, droplet 
size, equivalence ratio and amount of prevaporized fuel. These predictions 
were made with two existing spark ignition models which used a 
characteristic time approach. An ignition frequency model was developed 
which extended the characteristic time models to account for the 
probabalistic nature of the ignition process. 

Chapter 2 includes background information concerning the ignition of 
combustible mixtures, as well as a review of some previous work in the 
area. The spray characterization and ignition facilities and experimental 
results are described in Chapters 3 and 4, respectively. Analytical 
calculations and predictions are presented in Chapter 5. The conclusions of 
the study and recommendations for future work comprise Chapter 6. 


CHAPTER 2 


BACKGROUND AND LITERATURE REVIEW 

An electric spark ignites a combustible mixture by means of heat addition 
followed by flame propagation. Therefore, insight into the spark ignition 
process can be gained by studying the related processes of thermal ignition 
and flame propagation. Experimental and theoretical studies done in these 
areas will be reviewed in the first two sections of this chapter. 

Experimental spark ignition studies have been performed with both 
homogeneous and heterogeneous mixtures. Those studies will be reviewed in 
the third section of this chapter. The discussion will begin with a brief 
overview of the homogeneous studies coupled with a review of spark 
characteristics and their effects on ignition. This will be followed by a more 
detailed review of the heterogeneous ignition studies. 

Theoretical work on the spark ignition of sprays will be covered in the 
final section of this chapter. Global ignition models will be reviewed, as well 
as attempts to model the local formation and growth of the spark kernel and 
its effect on spray droplets within and surrounding it. 

2.1 Thermal Ignition Studies 

The ignition characteristics of single isolated fuel droplets in a hot 
environment was studied experimentally by Saitoh st li- (1982). They 
measured the ignition delay time as a function of droplet size and air 
temperature for n-heptane and n-hexadecane droplets ranging from 0.5-2. 5 



mm. A minimum droplet size tor ignition was observed, and below this 
minimum the droplet could completely vaporize without igniting. This 
limiting drop size, as well as the ignition characteristics of fuel droplets in 
general, depended strongly on the fuel type and ambient conditions. 

Law and Chung (1978, 1980) theoretically studied the thermal ignition of 
single droplets in a hot ambient environment. They found that the ignition 
delay time could be represented in terms of an ignition Damkohler number. 

The thermal ignition of monosized droplet streams in a hot environment 
was Investigated by Sangiovanni and Kesten ( 1977) . They studied the effect of 
droplet interactions, as controlled by the droplet size (200-300 Mm) and 
interdroplet spacing (2-200 droplet diameters), on the ignition delay time of 
furfuryl alcohol droplets. Theoretically, they determined ignition delay 
times for the limiting cases of an isolated droplet and a cylindrical filament 
of fuel. Results showed that the ignition delay time increased substantially 
as the interdroplet spacing was decreased. This droplet interaction effect 
increased with decreasing droplet size, decreasing ambient temperature and 
increasing fuel volatility. 

Aggarwal and Sirignano (1985) developed a numerical model for the 
ignition of idealized, quiescent monosized droplet arrays by a heated wall. 
They predicted optimum values of droplet size and overall equivalence ratio 
for ignition. The optimum droplet size depended on the overall equivalence 
ratio and was attributed to the distribution of fuel vapor in the ignition zone, 
which was affected by the total droplet surface area. Similarly, the optimum 
equivalence ratio was a function of the initial droplet size and fuel volatility. 
The authors also showed that ignition of droplet arrays (and fuel sprays) was 
a probabalistic phenomenon, since it depended on the location of the droplets 
relative to the heat source, which was random. They subsequently extended 
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this work to polydisperse sprays (Aggarwal and Sirignano, 1986) and 
concluded that the area mean diameter better characterized polydisperse 
sprays than the Sauter mean diameter. Aggarwal and Nguyen (1987) 
performed a numerical simulation of the ignition of monosized droplet arrays 
flowing in a tube. Their ignition source was a localized heat source (flame 
kernel), and they determined minimum ignition energies and ignition delay 
t times for n-hexane and n-decane droplets ranging from 30-150 urn. Results 

' showed that the minimum ignition energy increased with increasing droplet 

i 

size, increasing mixture velocity and decreasing fuel volatility. They also 
indicated an optimum droplet size (minimum ignition delay time) for 

i 

i n-hexane, and that this optimum size increased with increasing mixture 

| velocity. 

2.2 Flame Propagation and Transition Region Effects 

| Burgoyne and Cohen (1954) studied the effect of droplet size on flame 

propagation through monodisperse tetralin aerosols. They observed that 

premixed, prevaporized combustion dominated for sprays with droplet sizes 

% 

below 10 urn. Above 30 jam the droplets burned individually in their own 
diffusion flame envelope, with burning droplets igniting adjacent droplets. 
For the combustion of sprays with intermediate droplet sizes (10-30 pm), 
individual droplets encased in flame envelopes were observed superimposed 
on the gaseous laminar flame structure. Sprays which exhibited this 
transitional combustion behavior, both heterogeneous and homogeneous, 
were said to be in the ’transition region’. They also observed that the 
presence of droplets extended the lean prevaporized flammability limits to 
leaner equivalence ratios. They explained this with the observation that even 
at lean overall equivalence ratios, fuel droplets supported a diffusion flame 



burning In stoichiometric proportions. Conversely, a lean prevaporized 
premixed mixture would never have regions of stoichiometric mixture 
strength. 

Polymeropoulos and Das (1975) observed a maximum flame speed for 
kerosine sprays at about 30 Mm SMD . This was further evidence that flame 
propagation was enhanced by transition region effects. Hayashi liii- ( 1976, 
1981) studied flame speeds in a closed combustion bomb, which was subject 
to large pressure increases as the reaction progressed. However, they 
found that the flame speed increased as drop size increased from 16-40 urn, 
provided that the mass of liquid fuel droplets was greater than the mass of 
fuel vapor present. They also showed that the presence of fuel droplets 
could increase the flame speed compared to that for the fully vaporized 
case. 

Mizutani and Nakajima (1973) observed that the addition of kerosine 
droplets to a gaseous propane-air mixture increased the flame speed at lean 
overall equivalence ratios. The optimum drop size for enhancement was 
found to be 40 um. They explained this behavior by proposing that the 
droplets wrinkled and lacerated the flame surface, thereby increasing the 
flame speed. They also suggested that the drops served as high temperature 
ignition sources, extending the flammability limits and accelerating the 
burning velocity of adjacent flame elements. The burning velocity was 
increased because the flame propagated through regions of optimum fuel to 
air ratio formed around the evaporating drops. In addition, the burning 
droplets caused local gas expansions, generating turbulence which further 
intensified the transport processes. By repeating their experiments with 
vaporized kerosine, they confirmed that it was the two-phase nature of the 
flow that was responsible for the observed behavior. 


In addition to enhanced flame speeds and broadened flammability limits, 
several other combustion phenomena have been observed for sprays in the 
transition region. Nizami sill. (1982) and Sarv (1985) measured minima in 
NO x emissions for sprays in the 40-60 Jim size range. Chan (1982) and Singh 
(1986) both observed optimum drop sizes for spark ignition. These ignition 
studies will be reviewed in detail below. 

2.3 Experimental Spark lonition Studies 

2.3.1 Soark Chara cteristics a nd Homogene ous Spark lonition 
Some of the earliest experimental ignition work was performed by Lewis 
and Von Elbe (1961) and Litchfield (1960), who independently determined the 
minimum spark energy required for ignition of mixtures of 8.5% methane in 
air to be 0.28 mJ. In both studies, purely capacitive sparks with durations 
of less than 1 us were used as the ignition source. The spark energy (E sp ) 
was assumed to be equal to the stored electrical energy of the capacitor, 
given by: 

E sp = 0 .5CV 2 (2-1) 

where C is capacitance and V Is voltage. 

Rose and Priede (1958) investigated the minimum ignition energy of 
hydrogen-air mixtures using a capacitive discharge ignition circuit with 
series resistance (RLC circuit) to produce their ignition sparks. They found 
that by increasing the resistance in their circuit while holding the spark 
energy constant, the energy required for ignition decreased. Increasing the 
circuit resistance increased the duration of the spark. This decreased the 
strength of the shock wave produced by the spark, thereby increasing the 
proportion of the spark energy available for heat addition to the mixtures. 



Since some of the stored capacitor energy was dissipated in the circuit 
resistance, spark energies had to be calculated by integrating the product of 
the voltage (V) and current (I) across the spark gap over the duration of the 
spark: 


E S p*oJ tvldt <2-2 > 

They also noted that changing the electrode spacing (spark gap width) 
affected their ignition energy results and that this parameter should be 
optimized for each condition tested. 

In a detailed investigation of spark discharge characteristics Maly and 
Vogel (1978) studied the three modes of capacitive discharge: breakdown, 
arc and glow. Breakdown discharge occurs when there is very little or no 
resistance in series with the capacitor, resulting in spark durations that are 
very short, typically less than I ms. Arc and glow discharges have longer 
durations (10-300 ms and greater than 300 ms, respectively) due to 
increased series resistance. Maly and Vogel concluded that breakdown was 
the most efficient discharge mode because very little energy loss due to 
conduction of heat from the spark to the electrodes occurred, while 
substantial conduction losses occurred during both arc and glow discharges. 
The energy losses due to radiation were also found to be negligible for all 
three modes. These findings confirmed earlier studies on spark discharges 
by Roth it M- (1951). 

In a study similar to that of Maly and Vogel, Kono at $1. , (1977) 
concluded that arc discharge was the most efficient mode of energy transfer 
for the spark ignition of quiescent propane-air mixtures. They observed that 
under the condition of optimum spark duration and gap width, arc discharges 
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resulted in the lowest ignition energy of the three modes. 

While investigating the spark ignition of flowing propane-air mixtures, 
Swett (1956) found the optimum spark duration to be about 100 ms. Ballal and 
Lefebvre (1975) also observed optimum spark durations of 60-100 ms for 
flowing propane-air mixtures, as did Ziegler (1984) for methane-air 
mixtures. In addition, Ballal and Lefebvre (1975) found that the minimum 
ignition energy decreased with decreasing electrode diameter due to 
decreasing electrode surface area available for conduction losses. 

In summary, arc discharge has been found to be the most efficient form 
of energy transfer for the spark ignition of combustible mixtures, provided 
the discharge occurs under optimum conditions of spark duration and spark 
gap width to minimize energy losses. In reality, all arc discharges are 
initiated by a breakdown discharge, needed to ionize the medium in the spark 
gap and provide a conducting path for the arc. However this breakdown phase 
accounts for a negligible amount of the spark duration and total spark energy 
relative to the arc discharge. 

2.3.2 Heterogeneous Spark Ignition 

Spark ignition of polydisperse fuel sprays has been studied extensively by 
Lefebvre and co-workers. Rao and Lefebvre (1976) measured the minimum 
ignition energy (E m j n ) of polydisperse kerosine sprays with Sauter mean 
diameters (SMD) ranging from 30-100 Mm and equivalence ratios (0) ranging 
from 0. 4-1.0. They observed that droplet size had the strongest effect on 
ignition, with E m j n Increasing steadily with increasing SMD. They also found 
that, for a fixed droplet diameter, E min decreased steadily with increasing 
0 . 

Continuing this work, Ballal and Lefebvre (1978) measured E m j n over 



similar ranges of SMD and 0 for polydisperse sprays of six different fuels 
ranging in volatility from iso-octane to heavy fuel oil. The observed effects 
of SMD and 0 on E m j n were similar to those noted above for kerosine. 
Moreover, the authors quantified these effects, showing that E mjn was 
proportional to SMD 3 and 0" 1 - 5 . In addition, E mtn was seen to decrease with 
increasing fuel volatility. They concluded that the most critical factor in 
determining ignition was the mass concentration of fuel vapor generated by 
the spark In the ignition zone. In the same study, Ballal and Lefebvre 
investigated the effect of pressure, mean velocity and turbulence intensity on 
ignition and found that E m j n increased with decreasing pressure, increasing 
flow velocity and increasing turbulence intensity. 

The first spark ignition study performed with monosized sprays was by 
Chan (1982), who studied the ignition of monodisperse tetralin sprays with 
drop sizes ranging from 8-32 pm and equivalence ratios of 0.4-I.0. He 
observed an optimum droplet size for ignition of about 15 jam and reasoned 
that as the drop size decreased from 30 to 15 urn the enhanced evaporation 
made ignition easier. As the drop size decreased below 15 pm, however, the 
spray approached a homogeneous quality and lost the benefit of individual 
droplet combustion, thus increasing the energy needed for ignition of lean 
mixtures. 

Singh (1986) performed ignition studies of monodisperse tetralin sprays 
and observed an optimum drop size for ignition of 22-26 urn, depending on 
stoichiometric conditions. He concluded that the amount of fuel vapor 
generated by the spark depended on the drop size. Sprays with drop sizes 
smaller than the optimum produced fuel vapor concentrations during the 
spark discharge that were richer than optimum for ignition, while drop sizes 
larger than optimum produced fuel vapor concentrations leaner than optimum 


for ignition. He also reported ignition results in terms of ignition frequency 
rather than minimum ignition energy, since there was no sharp delineation 
between spark energies which produced ignition and those which did not. 

2.4 Theoretical Spark Ignition Studies 

The analysis and modeling of spray formation, evaporation and 
combustion has been quite extensive. Two review papers which cover this 
area quite thoroughly have been published by Sirignano (1983) and Faeth 
(1979). 

Modeling of the spark ignition process itself was first performed by Lewis 
and von Elbe ( 1961 ) who related ignition energy to the laminar flame speed of 
the mixture. They assumed that the entire spark energy was instantaneously 
converted to thermal energy within a spark kernel. For successful ignition, 
this spark kernel should achieve some critical volume, or quenching 
distance, resulting in flame propagation. This critical volume was achieved 
when the heat release rate within the kernel equalled the rate of heat loss 
from it. 

More recently, Peters and Mellor (1980) used a characteristic time 
approach to the model ignition energy of quiescent sprays. They 
subsequently extended the model for use in gas turbine engines (Peters, 
1981; Peters and Mellor, 1982), predicting ignition energies as a function of 
equivalence ratio, SMD and pressure for quiescent and turbulent sprays. 
Their criterion for ignition was that the residence time of the fuel-air mixture 
at the spark gap must be longer than the time required for evaporation and 
chemical reaction. The residence time was the ratio of the quenching 
diameter to the mean velocity, and the evaporation and ignition delay times 
were given by the law and an Arrhenius expression, respectively. These 


models predicted ignition energies which correlated well with the 
experimental data of Ballal and Lefebvre (1977). Ball al and Lefebvre ( 1981 ) 
developed a similar model to predict Ignition energies for homogeneous and 
heterogeneous mixtures under flow conditions ranging from quiescent to 
turbulent. They related the time required for evaporation and chemical 
reaction to the time required for the hot kernel to be quenched by heat loss to 
the surrounding spray. Their ignition delay time was related to the flame 
speed of the mixture, rather than the activation energy used by Peters and 
Mellor. 

While the spray ignition models noted above have taken a global approach 
to spark ignition, other investigators have used a more localized approach 
by taking into account the time dependent nature of the spark discharge. 
Maly (1981) developed a temporally resolved thermal ignition model which 
took into account the unsteady nature of the ignition process in a detailed 
analysis of the different phases of spark discharge and ignition. He proposed 
that the energy available for ignition was concentrated in a narrow shell 
around the spark plasma. His criterion for ignition was that the expansion 
velocity of this outer shell must be equal to that of a self supporting flame 
front of the same size. 

Adelman (1981) developed a theory which related the time dependent 
energy input from a spark to the kernel expansion. He took into account the 
strong shock wave spherically expanding from the spark, followed by the hot 
expanding spark kernel. His criterion for ignition was that the spark kernel 
must grow to a critical size before its expansion velocity falls below a 
critical velocity. Using a similar approach, Singh (1986) applied a model of 
spark kernel growth, including shock wave effects, to his experimental 
ignition data for qualitative correlation. For sprays of very small droplets 
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(less than 10 urn), the model predicted that the local flow field generated by 
the spark discharge resulted in a droplet free annular region around the 
spark kernel. This fuel free zone acted as a barrier to flame propagation, 
explaining his experimental observation that sprays with very small droplets 
(6.7 jam) were very difficult to Ignite. 

Bradley and Lung (1987) developed a hydrodynamic model to predict the 
temporal development of the radial profiles of density, pressure, velocity 
and temperature from a spark discharge channel. They showed that the slow 
spread of thermal energy from the relatively narrow discharge channel was a 
limiting factor in spark ignition This spread could be accelerated by higher 
energy input in the early stages of the spark discharge, resulting in stronger 
shock waves, thermal waves and outward convection of energy. 


CHAPTER 3 


EXPERIMENTAL FACILITIES AND METHODS 

The basic spray generation and characterization facilities are described 
in this chapter, along with the spray ignition system and experimental 
procedures. Representative spray characterization results are presented 
and discussed in this chapter, while comprehensive spray characterization 
results appear in Appendix A. Actual spray operating conditions and ignition 
results will be reported and discussed in Chapter 4. 

3. 1 Sprav Generation 

3.1.1 Monodisperse Sprav Generation System 

A Berglund-Liu Vibrating Orifice Monodisperse Aerosol Generator, Model 
3050 (Berglund and Liu, 1973), was used to produce sprays of monosized 
droplets for this study. Figure 3.1 shows a photograph of a typical spray 
enlarged 20 times. The picture was taken as the spray exited the test section 
10 cm downstream of the point of generation. The curved white line is the 
bottom rim of the test section. (For reference, the thickness of the white 
line is about 0.5 mm. ) As seen, the droplets are monosized, well dispersed 
and randomly positioned. This generator has been successfully used 
previously for a number of monodisperse spray combustion studies (Nizami 
iLli. , 1978-1982; Sarv, 1985). The generator was mounted on a traversing 
mechanism which allowed manual translation of the unit in the x-y-z 
coordinates. 
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Figure 3. 1 Photograph of a Monodlsperse N-Heptane Spray 
Enlarged 20X (D 0 = 50 m; 0 O = 0.55; z = 10 cm) 


A cross-section view of the aerosol generator Is shown In Figure 3.2. 
Monosized droplets were generated by applying periodic disturbances from a 
function generator to a piezoelectric ceramic, which in turn exerted 
mechanical vibrations on a liquid fuel jet emerging from an orifice plate 
seated in the ceramic. The perturbed jet then broke into discrete droplets 
with a standard deviation in size of approximately 1% of the mean diameter. 
One droplet was created per cycle of disturbance, so the initial droplet 
diameter (D 0 ) could be calculated from the volumetric fuel flow rate ( Qf ) 
and the frequency of disturbance (f<j) and is given by: 

D 0 = (6Q f /rrf d ) 1/3 (3-1) 

The droplet diameter was controlled by varying the fuel flow rate and the 
frequency of disturbance. While the droplet diameter was not a direct 
function of the orifice diameter, a given sized orifice was only capable of 
producing monosized droplets within a relatively narrow 'monodisperse' 
range. Therefore, different diameter orifices were used to produce a broad 
range of monosized droplet diameters. Table 3.1 shows the monodisperse 
ranges obtainable for the different orifices and fuel flow rates used in this 
study. 

3.1.2 Fuel and Air Flow Systems 

Figure 3.3 shows a simplified schematic diagram of the interior of the 
spray generator. Also included are the radial and axial coordinate axes (r 
and z) used in this study. An infusion syringe pump (Harvard Apparatus 
Model 901) supplied liquid fuel to the spray generator. As the stream of 
uniform fuel droplets exited the orifice, it entered a turbulent jet of 
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Figure 3.2 


Cross Section View of the Monodlsperse 
Spray Generator 




ORIFICE 

FUEL 

FUEL 

MONODISPERSE 

DIAMETER 


FLOW RATE 

RANGE 

(Jim) 


(cc/m) 

(um) 

7.5 

METHANOL 

0.075 

16-20 

12.5 

METHANOL 

0.075 

20-30 

17 

N-HEPTANE 

0.18 

30-43 

22 

N-HEPTANE 

0.19 

45-57 

27 

N-HEPTANE 

0.19 

54-70 

17 

METHANOL 

0.35 

30-45 

22 

METHANOL 

0.37 

45-58 

27 

METHANOL 

0.375 

50-67 


Table 3.1 


Monodisperse Operating Ranges for Various 
Orifices and Fuel Flow Rates 




dispersion air which kept the droplets from coagulating. This air jet was 
formed by the dispersion air flow exiting the dispersion cup through a 1 mm 
diameter orifice. The dispersion cup was mounted so the dispersion orifice 
was 1 mm directly below the fuel orifice. 

The dispersion jet with droplets in suspension continued downward 
through the reducing section, where a co-flow of dilution air was entrained 
into the spray. This reducing section had smooth rounded edges to minimize 
flow recirculation. The spray then entered the test section, which was a 1 .6 
cm l.D. Pyrex tube 18 cm in length, and exited to the atmosphere. 

Both the dispersion and dilution air flow rates were set by electronic flow 
controllers (Tylan Model FC-260), and monitored with mass flow meters 
(Hastings Model ALL-5K) . The overall stoichiometry (0 O ) of the sprays was 
set by the relative flow rates of fuel, dispersion air and dilution air. For a 
given drop size spray, the fuel and dispersion air flow rates were fixed, and 
the equivalence ratio was varied by changing the dilution air flow rate. 

3. 1 .5 Fuel Prevaporization System 

In order to perform ignition energy measurements on prevaporized, 
premixed fuel/air mixtures, a fuel prevaporizing apparatus was added to the 
experimental facility. A schematic of this system is given in Figure 3.4. 

Liquid fuel and air entered a spherical 500 cc pyrex vessel, heated to 
approximately 150 *C, where the liquid fuel evaporated and mixed with the 
air. The mixture then flowed through a heated secondary mixing section 
constructed of 1.2 cm 1.0. stainless steel tubing 40 cm in length. This 
section had 4 stainless steel screens spaced evenly along its length to induce 
turbulence, thus ensuring uniform mixing of the vaporized fuel and air. The 
mixture then flowed into the dispersion air inlet of the aerosol generator. 
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Figure 3.4 Schematic Diagram of the Fuel Prevaporizing 
and Premixing Apparatus 




Since no dilution air flow was needed, the test section was attached directly 
to the outlet of the dispersion cup for the prevaporized ignition experiments. 
The test section used was identical to the one used for the spray 
experiments, as described in the previous section. 

3.2 Number Density Measurement 

A major portion of the spray characterization procedure involved 
specifying the local droplet number density at the spark gap, which was 
directly related to the stoichiometry. Spray number densities were 
determined from measurements made with a laser Doppler velocimetry (LDV) 
system. These measurements included droplet velocity and rate, air 
velocity and LDV probe area, and are detailed in the following sections. 
Representative results from these measurements are also included below, 
while comprehensive results appear in Appendix A. 

3.2.1 LDV System 

Mean droplet and air velocities, turbulence intensities and droplet rates 
were measured using a dual beam LDV system in the backsc3tter mode. A 
schematic of the system is shown in Figure 3.5. The basic system 
components were the laser (Spectra-Physics Model 124B Laser, Model 255 
Power Supply), optics (TSI Model 900 Series), photomultiplier (TSI Model 
962) and signal processor (TSI Model 1984 input conditioner and TSI Model 
1985 timer). The signal from the processor was sent through an analog to 
digital converter and analyzed with a DEC LSI-11 microcomputor. The 
FORTRAN and MACRO data analysis programs are listed in Appendix B. A 
complete description of the physical principles involved in LDV measurements 
can be found in a paper by Stevenson ( 1976). 







3.2.2 Sorav Number Density and Stoichiometry 

A spray characterization procedure was performed to specify the 
relative amounts of fuel and air in the spark gap region during ignition 
testing. The procedure involved making LDV measurements at points on a 
horizontal rectangular grid at the axial location of the ignition electrodes. 
Due to the difficulty in attempting LDV measurements through a cylindrical 
test section, a shorter test section cut off 3 mm above the laser beams was 
used for this procedure. This allowed LDV measurements to be made through 
air, with no obstruction from the test section. The basic procedure was as 
follows. 

Droplet rate (N) and velocity (V^) measurements were made in a 2 mm 
grid pattern over the entire cross section of the spray. The droplet number 
density (n) was determined locally in the sprays as follows: 

n = N/V a A p (3-2) 

where V a is the air velocity and A p is the probe area. V a was measured as 
described in section 3.2.4. A p was measured as described in Section 3.2.3. 
The number density was then converted to the local equivalence ratio by: 

0 = n(p f nD 0 3 /p a )/3(f/a) m>st (3-3) 

where Pf and p a are the densities of the liquid fuel and air respectively, D 0 is 
the initial droplet diameter and (f/a)m,st * s the stoichiometric fuel to air 
mass flow ratio. 

The initial droplet diameter and not the local diameter was used in order 
to account for all of the fuel (liquid + vapor) present locally. It was 


assumed that the evaporated fuel stayed in the radial vicinity of the droplets 
from which it evaporated, so that the radial profiles of fuel vapor matched 
the droplet number density profiles in the sprays. This was based on a 
characteristic mixing length analysis for a turbulent free jet (Hinze, 1981) 
from which radial mixing lengths of 1-2 mm were calculated for the sprays. 
Since these mixing lengths were very close to the average interdroplet 
spacings of the sprays (0.6-1 .5 mm), this assumption was reasonable. 

Typical radial profiles of number density and local equivalence ratio are 
given in Figure 3.6 for 50 urn methanol and n-heptane sprays with overall 
equivalence ratios of 0. 51 . Notice that while the equivalence ratio profiles of 
these sprays were very similar, the methanol number densities were roughly 
twice those of the n-heptane spray. This was because methanol, with fewer 
carbon atoms and a fuel-bound oxygen atom, required about twice the fuel to 
air ratio of n-heptane to achieve the same equivalence ratio. Equivalence 
ratio profiles for all sprays tested are shown in Appendix A. 

The spark gap equivalence ratio (0g) was defined as the average value of 
the local equivalence ratio over the width of the electrode spacing. Values of 
0 g are given in the next chapter in Tables 4.2 and 4.3 for all sprays tested, 
along with specific values of the spark gap widths used. Spark gap widths 
ranged from 2-7 mm, increasing with decreasing equivalence ratio. Due to 
the increasing amount of dilution air entrained into the dispersion jet, the 
equivalence ratio profiles became flatter in shape with decreasing overall 
equivalence ratio. This meant that larger gap widths corresponded to flatter 
equivalence ratio profiles. Therefore, the local equivalence ratio was 
relatively constant over the width of the spark gap for most cases. 

Figure 3.7 shows spark gap equivalence ratio (based on a 2 mm gap 
width) plotted as a function of the overall equivalence ratio for n-heptane 
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Figure 3.6 Radial Profiles of Droplet Number Density and 
Equivalence Ratio for Methanol and N-Heptane 
Sprays (D 0 = 50 jim; 0 O = 0.51; z = 10 cm) 
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Figure 3.7 Variation of Spark Gap Equivalence Ratio with 

Overall Equivalence Ratio for N-Heptane Sprays 
(D 0 = 50 jam; z = 10 cm; Spark Gap = 2 mm) 



sprays with 50 um Initial droplet diameters. As seen, 0g was larger than 
and increased linearly with 0 O . These trends were observed for all sprays 
tested. The uncertainty in 0g was determined from triplicate measurements 
of the same spray and found to be ± 5%. The symbol 0 will be used 
interchangeably with 0g subsequently, while the overall equivalence ratio 
will be denoted exclusively by 0 O . 

At this point, a slight digression is in order to report an interesting 
phenomenon observed during the spray characterization procedures. As will 
be discussed more fully in Section 3.2.4, the droplet velocities were greater 
than the air velocities for most of the spray cases studied. For these cases, 
there existed an instantaneous droplet number density (n inst ) which differed 
from the previously defined droplet number density. The instantaneous 
number density was a measure of the number of droplets per volume of space 
(not air), it can be thought of as a snapshot (frozen in time) of the droplets 
in space. It was determined using an equation similar to that used for the 
previous number density, but based on the droplet velocity instead of the air 
velocity: 


"inst sN/v d A p <3-4) 

where N is the LOV droplet rate, Vg is the mean droplet velocity and A p is the 
LDV probe area. 

A radial profile of the instantaneous droplet number density compared to 
the droplet number density is shown in Figure 3.8 for a 50 um initial droplet 
diameter methanol spray with an overall equivalence ratio of 0.42. As seen 
in Figure 3.8, the instantaneous number densities are less than the spray 
number densities. In the time frame of an ignition spark (*100 jjs), the 
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Figure 3.8 Radial Profiles of Instantaneous Droplet Number 
Density and Droplet Number Density for a Methanol 
Spray (D 0 = 50 ym; 0 O = 0.42; z = 10 cm) 





instantaneous number density is the more relevant quantity of the two to 
consider when characterizing the spray. However, for ignition to occur, the 
spark kernel must propagate into the surrounding mixture. The time frame 
of this propagation, based on the flame speed and width of the test section, is 
on the order of several milliseconds. During this time, the flux of both fuel 
and air across the kernel boundary must be considered. Therefore, the 
droplet number density based on the air velocity (n) is more relevant to 
ignition and will be used exclusively to characterize the spray equivalence 
ratios in the subsequent reporting and analysis of ignition data. 

3.2.3 LDV Probe Area 

The LDV probe volume is the ellipsoidal region where the two laser beams 
cross, creating the fringe pattern which produces Doppler bursts when 
traversed by a droplet. These Doppler bursts are analyzed to determine the 
droplet velocity, and are counted to determine the droplet rate or frequency. 
The LDV probe area (A p ) is the projection of the probe volume normal to the 
spray direction. The size of the probe area is a function of: the diameter of 
the laser beams; the diameter of the droplets traversing it; the index of 
refraction of the droplets; the laser wavelength and power; the LDV optics; 
and the LDV processor electronics. A p was determined in the following 
manner and subsequently used in the spray characterization procedure 
described in section 3.2.2. 

Measurements of droplet rate (N) were made in a grid pattern over the 
entire cross section of the spray, as seen in Figure 3.9. The droplet rate 
measured in the center of each 4 mm^ grid sector was assumed to represent 
the average droplet rate through the entire sector. The ratio of A p to the 
sector area (Aj) was equal to the ratio of measured drop rate (N pjj ) to the 
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Figure 3.9 Schematic Diagram of Grid Measurement Pattern 
for Spray Characterization and Probe Area 
Calibration Procedures 
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total droplet rate (Nj) traversing Aj. From conservation of droplets, the 
sum of Nj over the entire spray was equal to the droplet generation rate ( f^) . 
Therefore, the probe area was calculated by: 

A p = (ZjNpj)^ (3-5) 

These values were then used in Equations 3-2 and 3-4 to calculate the 
spray number densities. Measured values of the probe area are shown as a 
function of droplet diameter and fuel type in Figure 3. 10. As seen, the probe 
area varied from about 0.03 to 0.23 mm^, increasing with increasing droplet 
diameter. This was expected, since the intensity of scattered light increases 
with the diameter of the scattering particle squared. 

3.2.4 Air and Droolet Velocity 

For dilute sprays, such as those used in this study, it can be assumed 
that the fuel droplets have a negligible effect on the air flow (Faeth, 1979). 
Therefore, axial and radial air velocity profiles were measured initially over 
the entire range of air flow conditions and applied to all spray conditions 
tested in this study. 

Mean air velocities were measured with the LDV using the smallest 
droplets attainable (17 urn initial droplet diameter) as the flow seed. Some 
justification for this is required, as a 17 urn particle may seem too large to 
follow the air flow. Figure 3.11 shows the effect of droplet size on droplet 
velocity over the range of air flow rates used in the study. These data were 
taken for n-heptane sprays at the centerline of the test section at an axial 
distance 10 cm downstream of the dispersion orifice. As seen, the droplet 
velocities asymptotically approached a constant value as the initial droplet 


diameter decreased below 24 jim. This indicates that the smaller droplets 
have relaxed to the air velocity and can be used with the LDV to measure air 
velocities. It also indicates that the larger droplets (above 24 Mm) are 
travelling faster than the air at this location in the test section. To further 
illustrate the difference in velocity between the droplets and the air in the 
sprays, Figure 3.12 compares axial velocity profiles of 17 and 70 urn initial 
diameter droplets as they proceed downstream from the dispersion orifice 
along the test section centerline (recall that the 17 pm droplets are 
following the air velocity). Because the 70 pm droplets are exiting from a 
larger diameter fuel orifice, they start out more slowly than the smaller 
droplets. Further downstream however, the larger droplets maintain a 
higher velocity due to their inertia, while the velocity of the smaller droplets 
and air relaxes more quickly. Therefore, by the time the spray has reached 
the ignition electrodes, the larger droplets are traveling faster than the air. 
Droplet Reynolds numbers in the sprays, based on the droplet diameter and 
slip velocity, ranged from 0. 1-1.0, increasing with increasing droplet 
diameter. 

Typical profiles of the radial variation of mean air and droplet velocities 
across the test section are shown in Figure 3. 13. The droplet profiles were 
measured 10 cm axially downstream for methanol sprays with an overall 
equivalence ratio of 0.7 and initial droplet diameters of 35, 50 and 65 pm. 
Notice that both the air and droplet profiles are relatively constant across 
the center of the test section and tail off towards the edge. Again, the 
increase in droplet velocity with increasing size is clearly seen. It should be 
noted that the methanol sprays generally had higher droplet velocities than 
the n-heptane sprays. This is because methanol required roughly double the 
fuel flow rate of n-heptane to achieve the same equivalence ratio, as 
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Figure 3. 13 Radial Profiles of Mean Droplet and Air Velocity 
for Methanol Sprays of Varying Initial Droplet 
Diameter (0 O = 0.7; z = 10 cm) 



discussed previously. These higher fuel flow rates meant higher orifice exit 
velocities for the methanol droplets, which resulted in higher velocities 
downstream. It should also be noted that conservation of mass was obtained 
upon integrating the air velocity profiles across the section for all cases 
tested. 

RMS velocities were measured using the LDV system and converted to 
turbulent intensity by dividing by the mean velocity. Values of turbulent 
intensity for the droplet velocities ranged from 20-25 % at the center of the 
sprays (r = 0-5 mm) and increased slightly towards the edges. In general, 
the turbulent intensity of the droplet velocities increased with decreasing 
velocity and decreasing droplet size. 

3.3 Droplet Size Measurement 

3.3.1 Droplet Sizing System 

A Fraunhofer diffraction technique, based on a method developed by 
Dobbins et aL (1963) for polydisperse sprays, was used to measure the 
diameter of the droplets at the spark gap. A schematic of the droplet sizing 
system is given in Figure 3.14, and the method is described as follows. A 
collimated He-Ne laser beam was passed through the sprays, and the light 
diffracted by the droplets was focused by a 200 mm focal length collecting 
lens onto a photographic plate. If the diffracting droplets were monosized, 
the diffracted light appeared as a series of alternating light and dark 
concentric rings known as an Airy or Fraunhofer diffraction pattern (Dobbins 
£ili. , 1963). A photograph of the diffraction pattern produced by a 50 pm 
methanol spray is shown in Figure 3. 15. The bright center of the pattern is 
due to the incident laser beam. Four sets of dark and light diffraction rings 
are clearly visible around the center. The intensity of the bright rings 
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Figure 3. 14 Schematic Diagram of the Droplet Sizing 
Apparatus 





Figure 3. 15 Photograph of a Fraunhofer Diffraction Pattern 
produced by a 50 jim Methanol Spray 
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decreased radially outward, so there were actually several sets of larger 
rings not observed due to the resolution of the film. The size of the spray 
| droplets was obtained by measuring the radii of these light and dark rings, 

based on the following analysis. 

I The radial intensity distribution of a Fraunhofer diffraction pattern in the 

1 focal plane of the collecting lens is given by (Dobbins elai- * 1963): 

l(0)/l(O) = (2J,(a0)/a0) 2 ' (3-6) 

where 1(9 )/l (0) is the ratio of the intensity of diffracted light at angle 0 to 
the incident light intensity. Jj is a Bessel function of the first kind, a is the 
particle size number given by: a * ttD/X, where D is the diameter of the 
i diffracting droplets and X is the light wavelength (0.6328 urn). The 

^ collecting lens transforms the angular diffraction information to a radial 

coordinate system (with the small angle approximation) by the following 

' relation: sin0 * 0 * r/f, where r is the radius on the focal plane of the lens of 

! 

I the light diffracted at an angle 0, and f is the focal length of the lens. This 

allows a0 to be expressed as: 

a0 = rrDr/fX (3-7) 

The dark rings in the diffraction patterns occur when the diffracted light 
is zero, or at zeros of the Bessel function Jj (a0 ). These zeros occur at a0 = 

3.83, 7.02, 10.17, 13.32, 16.4, etc. The diameter of the diffracting 
droplets was directly determined by measuring the radius of the center of the 
dark rings and using the values of the Bessel function zeros above with 
Equation 3-7. The resulting equation for the droplet diameter is: 


(3-8) 


0 = BjXf/nrj 

where B, is the i ‘th zero of the Bessel function and r s is the radius of the 
center of the i ‘th dark ring. It can also be shown that the bright rings 
correspond to the zeros of a second Bessel function ^(oeO ) , occuring at ot9 = 
5.14,8.42, 11.62, 14.8, 17.96, etc. Therefore, the diameter of the spray 
droplets could also be determined with Equation 3-8 by measuring the radius 
of the light rings and using the zeros of 

The total number of rings measured for any given diffraction pattern 
varied from 4 to 8, increasing with increasing droplet diameter. The 
uncertainty in this droplet sizing technique was ± 1%, based on the variation 
in droplet diameter from multiple ring measurements. 

Using this method, the droplet diameter at the spark gap (Dg) was 
measured for each of the sprays tested in this study. Representative 
photographs of the diffraction patterns obtained for the entire range of fuel 
and droplet size sprays tested are given in Appendix A. Comprehensive 
droplet sizing results are given in the next chapter in Tables 4.2 and 4.3. 
For a given fuel and droplet size spray, there was less than 5% variation in 
the measured spark gap droplet diameter over the entire range of 
equivalence ratios. 

Figure 3.16 shows the axial variation in droplet diameter for methanol 
and n-heptane sprays of 50 urn initial droplet diameter and overall 
equivalence ratio of 1.0. The n-heptane droplets evaporated faster due to 
their greater volatility. The n-heptane sprays also had lower droplet 
velocities than the methanol sprays, as discussed in the previous section. 
This meant longer residence times between generation and ignition, which 
also led to enhanced evaporation. 
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Figure 3. 16 


Axial Variation in Droplet Diameter for 50 urn 
Initial Droplet Diameter N-Heptane and Methanol 
Sprays (r ■ 0 mm) 



Droplet size measurements were performed at the centerline of the test 
section. Since the droplet sizing technique made a line of sight average 
through the entire spray, no effect of radial position was expected. This was 
confirmed by measuring the drop size of a 50 urn n-heptane spray at three 
different radial positions (-4, 0 and 4 mm). The measured values were all 
within 1 %. 

3.3.2 Extent of Prevaporization 

The amount of evaporated fuel present in a spray is a critical factor in 
ignition (Ballal and Lefebvre, 1976). Increasing the extent of 
prevaporization not only decreases the energy required to evaporate the fuel 
droplets in the spark kernel, it also increases the rate of propagation, or 
growth, of the kernel into the unburned mixture. Therefore, an accurate 
measure of the extent of prevaporization was necessary for the sprays tested 
in this study. 

The extent of prevaporization (0) was defined as the ratio of the mass of 
evaporated fuel to the mass of total fuel present. For the prevaporized, 
premixed cases, 0 was equal to one. For the sprays, 0 was determined at 
the spark gap based on the initial droplet diameter and the measured droplet 
diameter at the spark gap. Values of 0 are presented in the next chapter in 
Tables 4.2 and 4.3 for the sprays tested. 0 was kept relatively constant 
(within 5%) for each fuel. This was done by varying the axial location of the 
ignition electrodes. For example, the smaller droplet size sprays, which 
evaporated more quickly, were ignited at at smaller axial distance 
downstream. 
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3.4 Sprav Ignition 

3.4.1 Spark Generation System 

Ignition sparks wore delivered to the sprays by electrodes which 
protruded into the test section at axial distances 7.5, 10 and 12.5 cm below 
the dispersion orifice. The electrodes were placed with their axes normal to 
the spray direction, and were mounted on a micrometer traversing assembly 
which was used to adjust the spark gap width. The electrodes were 5 mm 
diameter stainless steel rods, with the spark producing ends machined to tips 
of 1 mm diameter and 6 mm length. The high voltage side electrode tip was 
pointed while the ground side tip was flat, as seen in Figure 3. 17. The tips 
were bent upward approximately 15* to prevent accumulation of fuel droplets 
on them. 

A capacitive discharge spark generation system produced the ignition 
sparks at the electrodes. The system was designed to provide independent 
control of the energy and duration of individual sparks, and consisted of a 
modified RC circuit as shown in Figure 3. 18. 

Sparks were generated in the following manner. The capacitor 
(0.5-50nF) was charged to 20 kV by closing the - charging switch for 2 
seconds. The charging switch was then opened. The spark was generated 
immediately afterward by closing the discharge switch, which caused the 
capacitor to discharge through the resistor (1-1000 kQ) and across the 
electrode gap to ground. The two switches were high voltage solenoid 
switches, and were controlled by timing relays which allowed repetition of 
the charge/discharge cycle. The time interval between sparks could be 
varied from 3-15 seconds. The charging resistor was used to limit the 
current during charging, while the shunt resistor provided an alternative 
path to ground in cases where the discharging current could not jump the 
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Figure 3. 18 Schematic Diagram of the Spark Generation Circuit 



spark gap. 

The spark energy and duration were determined by the values of the 
capacitor and the resistor in the circuit. To measure the spark energy and 
duration, a high voltage probe was connected across the electrodes, and an 
inductive current probe was attached to the ground-side electrode. The 
resulting voltage and current traces of the spark were recorded on a two 
channel digital storage oscilloscope. Voltage, current and energy traces for 
a typical spark are shown in Figure 3. 19. As seen, after the initial spike of 
the breakdown voltage, the spark voltage was relatively constant throughout 
the spark duration, while the spark current and energy decayed 
exponentially. The energy of a given spark was determined from: 

E = oP V(t) Kt) dt (3-9) 

where V is voltage, I is current, t is time. The spark duration, t, was taken 
as the time required for the spark current to decay to 5% of its initial value. 

A series of calibration measurements was performed to determine the 
spark energy levels of the system described above. This spark energy 
calibration procedure was done over a wide range of spark gap widths (2-8 
mm) and capacitance (0. 1-50 nF) for three discreet values of spark duration 
(30, 60 and 1 00 us). These durations were chosen to cover the range of 
optimum spark durations observed by previous researchers for similar spray 
conditions, including Rao and Lefebvre (1976), Ball a 1 and Lefebvre (1978) 
and Chan and Polymeropoulos (1982). Typical results of these calibration 
measurements are given in Figure 3.20, which shows spark energy as a 
function of spark gap width for the range of capacitors and a duration of 60 
us. Each data point is the result of two sets of 15 spark energy 
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measurements, taken on separate days. The standard deviation of each set 
of 15 measurements was less than 10% of the mean, while the day to day 
variation in the mean spark energy values was less than 5%. As seen in 
Figure 3.20, the spark energy varied logarithmically with the spark gap width 
for a given capacitance, which is seen as a linear relationship on the 
semi-log plots shown. This calibration procedure was performed in air, but 
it was subsequently observed that the composition of the mixture in the spark 
gap affected the measured spark energies by altering the resistance of the 
spark gap. Therefore, spark energy levels were measured for each 
individual ignition case studied, as will be described in the following section. 

3.4.2 Minimum Ignition Energy Measurement 

Spark ignition is a probabilistic phenomenon, as was discussed in Chapter 
2. No sharp boundary was observed between spark energies which ignited a 
given mixture and those which did not. Rather, a range of spark energies was 
seen for which increasing the spark energy resulted in an increasing 
probability or frequency of ignition. An example of this is given in Figure 
3.21, which shows the effect of spark energy on ignition frequency* for an 
n-heptane spray (0g=O.88, 0g=4 1 urn). Notice that the ignition 

frequencyrose sharply initially, then tailed off as it assymptotically 
approached an ignition frequency of 1.0 (100%). Similar behavior was 
observed for all of the spray and prevaporized ignition cases studied. 


* The term 'ignition frequency' was employed based on previous usage in 
the literature. It should be noted, however, that this quantity is not a true 
frequency, but represents the incidence or probability of ignition. 
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Each point in Figure 3.21 was the result of ignition observations and 
spark energy measurements for one set of 50 single sparks. The interval 
between sparks varied from 5 to 20 seconds, increasing with decreasing 
equivalence ratio. This interval was always sufficient to allow the 
temperature and spray conditions in the test section to return to normal after 
a successful ignition. The criterion for ignition was that a visible flame 
propagate at least 5 cm into the mixture. The spark energy was measured 
simultaneously with the ignition observations and averaged for 20 of the 50 
sparks, as described in the previous section. A minimum of two sets of 50 
sparks was done at each energy level. The values of spark energy and 
ignition frequency for both sets were within 7% of the mean for all cases 
studied. 

The minimum ignition energy (E m j n ) was defined as the spark energy 
level which produced an ignition frequency of 0.5 (50% ignition), as shown in 
Figure 3.21. The 50% level was chosen to account for the probabilistic 
effects of ignition, as will be described in Chapter 5. The uncertainty in E mjn 
was calculated from the 95% confidence intervals of repeated sets of spark 
energy and ignition frequency measurements and found to be ± 10%. 

3.4.3 Soark Parameter Optimization 

The ignition frequency curve shown in Figure 3.21 was obtained with fixed 
values of spark gap width and spark duration. As discussed in Chapter 2, 
both of these parameters affect ignition. If the spark gap width is too small, 
excessive heat is lost to the electrodes through conduction. If it is too 
large, the spark kernel becomes larger than necessary for ignition, 
resulting in excessive heat lost to the surrounding mixture through 
convection. If the spark duration is too short, excessive energy is lost in the 



accompanying shock, wave, whose strength increases with decreasing spark 
duration. If the duration is too long, spark energy is still being input after 
ignition has occurred, resulting in convection losses to the mixture. 
Therefore, optimization of these spark parameters involves minimizing 
energy losses and must be done to obtain the lowest possible minimum 
ignition energy for each Ignition case studied. 

The procedure used to optimize the spark duration and gap width for a 
given spray case was as follows. The spark energy level (capacitor) which 
produced ignition frequencies in the 40% to 60% range was determined. With 
this capacitor, one set of 50 sparks was performed at each of three durations 
(30, 60 and 100 us) over a range of gap widths. Starting at a small value, 
the gap width was increased in steps of 1 mm. For a given duration, the 
ignition frequency increased with increasing gap width, reached some 
maximum, then decreased. This was repeated for a larger (or smaller) 
capacitor in order to obtain ignition frequencies both above and below 50%. 
The optimum values of duration and gap width were then easily determined as 
those which resulted in the lowest 50% energy. For cases where the optimum 
duration was 30 or 100 us, durations of 15 and 200 us respectively were 
tested to confirm an optimum value. 



CHAPTER 4 


EXPERIMENTAL RESULTS AND DISCUSSION 

The effects of droplet size, equivalence ratio and fuel properties on the 
minimum ignition energy of n-heptane and methanol sprays are reported in 
this chapter. Ignition energy results for prevaporized n-heptane and 
methanol mixtures are also reported to represent the lower limit of droplet 
size. Valuable insight into the ignition phenomenon is gained by comparing 
the spray and prevaporized results as a function of both equivalence ratio and 
droplet diameter. 

4.1 Sprav Ignition 

4.1.1 Operating Conditions 

The generation parameters and flow rates used to produce the n-heptane 
and methanol sprays ignited in this study are listed in Table 4.1. Also 
included are initial droplet diameters and overall equivalence ratios. As 
seen, air flow rates for the methanol and n-heptane sprays were similar, but 
methanol sprays required roughly double the fuel flow rates to achieve the 
same range of stoichiometry because of the fuel bound oxygen. For a given 
fuel and drop size spray, the dispersion air flow rate was set just high 
enough to adequately disperse the droplets. For a given drop size spray, the 
equivalence ratio was varied by adjusting the dilution air flow rate. 
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Table 4.1 


Spray Generation Parameters and Operating 
Conditions 


Tables 4.2 and 4.3 contain values of the relevant spray parameters 
measured at the spark gap, as described in Chapter 3, for n-heptane and 
methanol sprays, respectively. Also included are the optimized ignition 
parameters and measured minimum ignition energies, which will be discussed 
in the next section. 

While the tables are straightforward, several points should be noted. 
The axial location (z) of the ignition electrodes was varied as a function of 
droplet size, with the smallest sprays being ignited at the smallest axial 
distance. This was done for several reasons. The smallest droplets (30 and 
33 wm) spread to a greater degree than the larger droplets as they 
progressed downstream. This resulted in droplet impaction, beginning about 
10 cm downstream, on the test section walls. To solve this problem, the 
ignition electrodes were moved upward to an axial location of 7.5 cm where 
the sprays just filled the test section with no impaction. Conversely, the 
largest sprays (53 and 57 urn) did not spread sufficiently by the time they 
had reached 10 cm downstream. For these larger droplets the electrodes 
were moved to 12.5 cm where the sprays were more uniform across the test 
section. This variation in electrode positioning also had the effect of 
minimizing the difference in the extent of prevaporization between the 
different sized sprays. This resulted in values of Q within 5% for a given 
fuel. 

All sprays were generated with liquid fuel and air at 20 *C. The 
temperatures listed in Tables 4.2 and 4.3 were measured at the spark gap 
with a type K thermocouple and represent a bulk spray temperature. 
N-heptane sprays temperatures ranged from 8 to 10 *C while methanol spray 
temperatures ranged from -1 to 2 *C. The methanol sprays had lower 
temperatures because of methanol’s higher latent heat of vaporization. 
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Table 4.2 


N-Heptane Spray Characterization and Ignition 
Results 
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Table 4.3 Methanol Spray Characterization and Ignition 
Results 


cn 
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Greater heat transfer from the surrounding air to the droplet was required 
for evaporation, lowering the spray temperature. Both the lower spray 
temperatures and the higher latent heat decreased the volatility of the 
methanol sprays, resulting in their lower extent of prevaporization compared 
to n-heptane sprays. 

Mean values of the spray number density at the spark gap varied from 
0.5-4. 1 mm~3 for n-heptane sprays and 1.0-13 mm - ^ for methanol sprays, 
increasing with increasing equivalence ratio and decreasing droplet size. 
Mean values of the interdroplet spacing at the spark gap varied from 20-32 
droplet diameters for n-heptane sprays and 15-25 droplet diameters for 
methanol sprays, increasing with decreasing equivalence ratio and 
decreasing droplet size. 

4.1.2 Spark Duration and Gao Width 

Minimum ignition energies were measured at optimum values of the spark 
duration and spark gap width for each spray tested, as described in Section 
3.5.3. These optimum values are listed in Tables 4.2 and 4.3. Optimum 
values of both the spark duration and gap width decreased with increasing 
equivalence ratio. Optimum durations ranged from 60 to 100 us for both 
n-heptane and methanol sprays. Optimum gap widths ranged from 2-7 mm 
for n-heptane sprays and 4-8 mm for methanol sprays. These values are in 
the same range found by previous researchers in similar studies (Rao and 
Lefebvre, 1976; Ballal and Lefebvre, 1978; Chan, 1982; and Singh, 1986). 

Varying the spark gap width was found to have a large effect on E m j n , 
while varying the spark duration had little effect. Figure 4.1 shows the 
effect of spark gap width and duration on the minimum ignition energy of a 41 
urn n-heptane spray with an equivalence ratio of 0.69. As seen, varying the 
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gap width from 2 to 6 mm had a pronounced effect on E m j n , with an optimum 
value of 4 mm clearly evident for each duration. Varying the spark duration 
from 30 to 100 us had a lesser effect on E m j n , but an optimum duration of 60 
US was evident. Similar behavior was observed for all of the sprays tested. 

4.1.3 Effect of Droplet Size and Equivalence Ratio 

The variation in ignition frequency with spark energy was determined as 
described previously for the sprays listed in Tables 4.2 and 4.3 . Figure 4.2 
shows typical ignition frequency curves for 56 um n-heptane sprays ranging 
in equivalence ratio from 0.59 to 1 .80. Notice that as 0 decreased from rich 
to lean, the curves became less vertical and began to tail off towards higher 
spark energies as they approached 100% ignition. Similar behavior was 
observed for all cases tested. Ignition frequency curves for all cases are 
shown in Appendix A. 

The minimum ignition energy (E m j n = 50% ignition) was determined from 
these ignition frequency curves as described in Section 3.5.2. The effect of 
equivalence ratio and droplet diameter on E m j n is shown in Figure 4.3 for 
n-heptane and methanol sprays. As can be seen, both- sets of curves are 
very similar qualitatively. For a given droplet size and fuel spray, E mm 
decreased sharply with increasing lean equivalence ratio, but began to level 
off as the sprays became richer. E m j n continually decreased with increasing 
equivalence ratio and no optimum value of equivalence ratio for ignition was 
observed for any of the sprays tested. The lean ignition limit was determined 
to be about 0 = 0.4 for both n-heptane and methanol sprays. For any given 
equivalence ratio and fuel spray, E m j n decreased with decreasing droplet 
size. 

Similar behavior of E m j n with respect to droplet size and equivalence 
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Figure 4.3 Effect of Equivalence Ratio and Droplet Diameter 
on the Minimum Ignition Energy of N-Heptane 
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ratio has been observed in previous ignition studies with sprays in the same 
size range (Rao and Lefebvre, 1976; Ballal and Lefebvre, 1978 and 1981 ) . The 
decrease in E m j n with decreasing droplet size and increasing equivalence 
ratio can be explained as follows. For a given equivalence ratio, the smaller 
sprays had more total surface area of fuel which decreased the time required 
for the droplets to vaporize and burn in the spark kernel. The larger 
droplets took longer to vaporize and burn in the spark gap, resulting in 
energy losses to the unburned mixture. Therefore, more ignition energy was 
required for larger droplet size sprays. For a given droplet size, increasing 
the equivalence ratio improved ignition by increasing the amount of fuel in and 
around the spark kernel. This effect became less pronounced as the sprays 
became richer. 

4.1.4 Effect of Fuel Properties 

As seen in Figures 4.3, the methanol sprays required roughly 3 times 
more spark energy for ignition than the n-heptane sprays. This can be 
explained by examining the relative volatilities of the two fuels. Previous 
researchers (Ballal and Lefebvre, 1978; Peters and Mellor, 1980) have 
shown that fuel volatility is a critical factor in spray ignition, with E m j n 
decreasing with increasing volatility. The volatility of the fuel in the spark 
kernel during ignition can be quantified in terms of a combustion mass 
transfer number (B) which may be defined as the ratio of the energy available 
for evaporation to the energy required for evaporation. Values of B (Kanury, 
1975) were 5.82 and 2.70 for n-heptane and methanol, respectively. 
N-heptane's higher volatility, more than twice methanol's, accounted in a 
large part for the better ignitability of n-heptane sprays. 

Increasing the extent of fuel prevaporization prior to ignition (Q) has 


also been shown to reduce ignition energy requirements (Ballal and Lefebvre, 
1981). Values of 0 for the methanol sprays were about 25% lower than for 
the n-heptane sprays, as seen in Tables 4.2 and 4.3. This difference in 0 
also accounted in some part for methanol's higher ignition energies. The 
relative importance of fuel prevaporization and volatility on ignition will be 
further examined analytically in Chapter 5. 

Several differences between n-heptane and methanol combustion were 
observed in the types of flames produced by ignition. For very lean cases (0 
< 0.5) the n-heptane flames were faint blue in color. As 0 increased beyond 
0.5, they became a brighter orange color. This indicated increasing soot 
formation and oxidation. In fact, the soot had to be cleaned from the 
electrodes often while testing n-heptane sprays. The methanol flames were 
a very faint blue color at lean equivalence ratios and became faint orange 
with increasing 0, indicating very little sooting. This was confirmed by the 
fact that the electrodes showed no soot accumulation during methanol 
ignition. 

4.2 Prevaporized Ignition 

4.2.1. Operating Conditions 

Ignition energy experiments were performed on prevaporized, premixed 
flowing n-heptane and methanol mixtures over the full range of equivalence 
ratios attainable. The fuel and air flow rates used to produce these mixtures 
are listed in Table 4.4, along with the measured minimum ignition energies 
and optimized spark parameters. The ignition electrodes were located 10 cm 
downstream axially for all prevaporized experiments. The mixture 

temperatures at the spark gap were 25 °C ± 1 for all cases. Equivalence 
ratios higher than 3.0 could not be tested for n-heptane mixtures because of 


MIXTURE PARAMETERS: 


IGNITION PARAMETERS 


N-HEPTANE 


0 
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OUR (us) 

GAP (mm) 

Emin (mJ) 
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0.19 

2750 
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9 

* 1000 

0.65 


2538 

60 

8 

20 

0.7 


2357 

30 

7 

2.9 

0.8 


2062 

30 

6 

0.70 

0.9 


1833 

30 

5 

0.40 

1.0 


1650 

30 

4 

0.34 

1.5 
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- 

- 

< 0.2 

2.0 


825 

- 

- 

< 0.2 
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30 

4 

0.43 

3.0 
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0 

Qf (cc/m) Qa (cc/m) 

DUR (us) 

GAP (mm) 

Emin (mJ) 

0.6 

0.38 3075 

60 

9 

* 1000 

0.65 

2798 

60 

8 

5.0 

0.7 

2559 

30 * 

5 

0.93 

0.75 

2352 

30 
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2171 

30 
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0.34 

0.85 

2012 

30 

2 

0.23 

0.9 

1870 

30 

2 

0.20 

1.0 

1628 

- 

- 
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Table 4.4 


Prevaporized Mixture Generation Parameters 
and Operating Conditions 


fuel condensation in the dispersion cup. For the same reason, equivalence 
ratios above 1.0 could not be tested for methanol mixtures. However, the 
range of prevaporized mixtures tested was sufficient to provide a lower limit 
(0 = 0) for the study of droplet size effects on ignition. 

4.2.2. Spark Duration and Gap Width 

The optimum values of spark duration and gap width are listed in Table 
4.4 for prevaporized n-heptane and methanol mixtures. With the exception 
of the very lean cases (0 = 0.6 and 0.65) the optimum duration was 30 us for 
all mixtures. These durations were shorter than the comparable spray 
cases, which was expected since no time was required for fuel evaporation 
during ignition. The optimum spark gap widths varied from 4-9 mm and 2-9 
mm for n-heptane and methanol mixtures, respectively. These values were 
comparable to those observed for spray ignition. The effect of varying spark 
duration and gap width on the minimum ignition energy of prevaporized 
mixtures was also similar to the behavior observed and reported for sprays. 

4.2.3. Effect of Equivalence Ratio 

The procedure for measuring the minimum ignition energy of the 
prevaporized mixtures was identical to the one used for sprays. The ignition 
frequency was determined as a function of spark energy and E m j n was the 
spark energy which produced 50% ignition. Ignition frequency curves for all 
cases tested are presented in Figures 4.4, 4.5 and 4.6 for lean n-heptane, 
rich n-heptane and lean methanol mixtures, respectively. The prevaporized 
ignition frequency curves were very similar in shape to those observed for 
sprays. 

For two of the lean cases (n-heptane and methanol, 0 = 0.6) E m j n could 
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only be estimated because the highest spark energies attainable with the 
ignition circuit, about 100 mJ, resulted in ignition frequencies less than 50%. 
Therefore only partial ignition frequency curves were obtained. E m j n was 
estimated by extrapolating these partial curves to the 50% line based on the 
shapes of the other curves. E m j n was estimated to be 1000 ± 300 mJ for both 
cases. For three of the richer cases (n-heptane, 0 * 1.5, 2.0 and 
methanol, 0 ■ 1.0) only an upper limit could be determined for E m j n because 
the lowest spark energies attainable, about 0.2 mJ, produced 100% ignition. 
E m in ’ s reported as < 0.2 mJ for these mixtures . 

The variation in minimum ignition energy with equivalence ratio is shown 
in Figure 4.7 for prevaporized, premixed n-heptane and methanol. An 
optimum equivalence ratio for n-heptane ignition, between 0 * 1.5 and 0 * 
2.0, is evident from this figure. This is in very good qualitative and 
quantitative agreement with data presented by Lewis and Von Elbe (1961), 
who reported an optimum equivalence ratio of about 1.8 corresponding to a 
minimum ignition energy of about 0.23 mJ for prevaporized n-heptane. Rich 
mixtures of methanol could not be tested due to the condensation problem 
mentioned previously, so no optimum equivalence ratio for methanol was 
observed. Prevaporized lean ignition limits were observed at equivalence 
ratios of 0.55 for n-heptane and 0.5 for methanol. 

4.2.4. Effect of Fuel Properties 

As seen in Figure 4.7, prevaporized methanol was easier to ignite than 
prevaporized n-heptane. This may seem to contradict the spray results, 
which showed that n-heptane sprays were more ignitable. However, the 
reason the methanol sprays were more difficult to ignite was because they 
needed more spark energy for fuel vaporization due to their lower volatility. 
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Volatility Is not a factor In the prevaporized case, and activation energy is 
the critical fuel property. Methanol has a lower activation energy than 
n-heptane, 41.3 cal/gmol as opposed to 60.5 cal/gmol (Kanury, 1975), 
which accounted for its better ignitability in the vapor state. 

4.3 Sprav and Prevaporized Comparison 

4.3.1. Extension of Lean Ignition Limits 

Figure 4.8 shows the spray and prevaporized ignition data plotted 
together for n-heptane and methanol. Extension of the prevaporized lean 
ignition limit to lower equivalence ratios was evident for both n-heptane and 
methanol sprays. This extension can be explained as follows. When the 
ignition spark occurs, temperatures in the spark path reach several thousand 
degrees Kelvin (Maly, 1978). Fuel droplets present in the immediate vicinity 
of the spark vaporize and burn, along with the prevaporized fuel, to form the 
spark kernel. Successful ignition occurs if the spark kernel is able to 
propagate into the unburned mixture surrounding it. In spray ignition, the 
kernel is surrounded by fuel droplets and a vaporized fuel/air mixture with a 
range of equivalence ratios. Even at lean overall qap equivalence ratios, 
localized regions that are richer and easier to ignite exist around the 
droplets. These regions of favorable stoichiometry accelerate the growth of 
the spark kernel into the surrounding mixture. Thus, the presence of fuel 
droplets enhances lean spray ignition. In contrast, prevaporized premixed 
mixtures have a uniform distribution of fuel vapor, so kernel growth is 
limited by the ignitability of the stoichiometry being tested. 

Similar extension of lean flammability limits has been observed in 
previous spray studies. For example, Burgoyne and Cohen (1954) found that 
the flammability limits were extended leanward for tetralin sprays with 
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droplets In the 10-40 urn range. Also, Mtzutant and Nakajlma (1973) obtained 
extension of the flammability limits by the addition of kerosine droplets to a 
propane air mixture. 

4.3.2. Optimum Droolet Size for Ignition 

The spray and prevaporized ignition energy data are given as a function of 
droplet size for various equivalence ratios in Figure 4.9 for n-heptane and 
methanol. Prevaporized data were plotted at zero diameter, the lower limit 
of drop size. The dashed lines are an extrapolation from the lowest drop size 
studied to the prevaporized case. The minima in the extrapolated region of 
the Ignition curves suggest the existence of an optimum droplet size for 
ignition of lean sprays, occuring below 30 um for both fuels. The existence 
of this optimum droplet size was indicated over a larger range of lean 
equivalence ratios for n-heptane (0.4-0. 8) than for methanol (0.4-0.65). 
The optimum diameter could not be quantified from this study, obviously. 

The existence of an optimum droplet diameter for ignition can be 
explained as follows. A successful spark ignition event can be considered as 
two separate processes. Initially, the spark has to vaporize and burn the 
fuel in its path to create a spark kernel. Secondly, this kernel has to 
propagate into the unburned mixture. The amount of energy needed to 
vaporize and burn a given amount of fuel in the spark gap decreases 
monotonically as the droplet size decreases to zero, due to the enhanced 
evaporative quality of the fuel. Therefore, the existence of an optimum drop 
size for ignition must be explained in terms of the second process, flame 
propagation. 

Flame propagation studies (Burgoyne and Cohen, 1954; Polymeropoulos 
and Das, 1976; Hayashi sial- , 1976; Mizutani and Nakajima, 1973) have shown 
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Figure 4.9 Effect of Droplet Diameter and Equivalence Ratio 
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that as drop sizes are decreased from large values to zero, sprays in the 
transition region (10-40 urn) achieve the maximum flame speed. This 
enhanced flame speed would also enhance spark kernel propagation rates, 
and ignition. Therefore, the optimum in drop size for ignition is probably 
caused by the increased flame propagation rates in the transition region, and 
should occur at the drop size where the flame speed is maximized. 

Previous researchers also have observed optimum droplet sizes for 
ignition. Chan (1982) and Singh (1986) found optimum sizes of 15 urn and 
22-26 urn respectively for ignition of tetralin aerosols at lean equivalence 
ratios, which agrees with the current findings. A numerical thermal ignition 
study by Aggarwal and Sirignano (1984) predicted an optimum droplet size 
for ignition, which increased with increasing equivalence ratio. This optimum 
was attributed to the distribution of fuel vapor in the ignition zone, which is 
affected by the total droplet surface area. 

4.4 Summary of Experimental Results 

To briefly summarize the major experimental results: 

1) Ignition was improved (required ignition energies were lowered) by 
decreasing the droplet size, increasing the equivalence ratio and increasing 
the fuel volatility of the sprays. 

2) An optimum equivalence ratio for ignition of prevaporized n-heptane 
was observed between 1.5 and 2.0; no optimum equivalence ratio was seen 
for spray ignition. 

3) An optimum droplet size for spray ignition, less than 30 urn, was 
suggested for lean sprays of both fuels. This optimum was related to the 
maximum flame speeds observed by previous researchers for sprays in this 


transition region size range. 

4) Extension of the lean prevaporized ignition limit was observed for all 
sprays tested. This was attributed to interdroplet regions of gas-phase 
stoichiometry more favorable to ignition created by the evaporating droplets. 


CHAPTER 5 


ANALYTICAL SPARK IGNITION STUDIES 

It has been shown that many factors affect the minimum Ignition energy of 
fuel sprays. Droplet size, equivalence ratio, extent of prevaporization and 
fuel properties are parameters that will be examined by two analytical 
models which are described in this Chapter. The models will be used to 
predict minimum ignition energies based on these and other parameters using 
a characteristic time approach. The results will give insight into the relative 
importance of these parameters on the spark ignition process, as well as the 
ability of characteristic time models to accurately predict ignition energies. 
One of the models will also be extended to simulate the randomness of the 
spray ignition process in order to explain the probabilistic nature of the 
experimental results. 

5. 1 Sprav Ignition Model (Ballal and Lefebvre) 

The ignition model developed by Ballal and Lefebvre (1981) takes a 
characteristic time approach, with the criterion for ignition being that the 
time required for the fuel to evaporate and chemically react in the spark 
kernel must be less than or equal to the time it takes the kernel to be 
quenched by thermal conduction and turbulent diffusion. A critical quenching 
diameter was calculated from this. The minimum ignition energy was the 
amount of energy required to raise this critically sized kernel to the 
stoichiometric adiabatic flame temperature. The fuel and air properties used 



in the model are listed in Table 5. 1 . 


5.1.1 Model Formulation 

The ignition model was developed to be used over a wide range of flow 
conditions. To correlate the present ignition data, the heterogeneous, 
quiescent and monodisperse case was used. (The model criterion for 
quiescent conditions was Req « 16. The droplet Reynolds numbers for the 
present sprays varied from 0.3-1 .2. ) 

The basic assumptions of the model were as follows. A spherical spark 
kernel was formed by the ignition spark. The temperature in this kernel was 
the stoichiometric, adiabatic flame temperature of the fuel. The air 
properties in the kernel were evaluated at a mean temperature of 1300 K. 
The basic premise of the model was given by: 

tq s t e + tc ( 5- 1 ) 

which equated the time required for droplet evaporation (t e ) and chemical 
reaction (t c ) to the quenching time (tq) of the spherical spark kernel. 

The quenching time was defined as the ratio of the heat capacity of the 
spark kernel divided by the average rate of heat loss from the kernel by 
conduction. This was given by: 

tq = c pa p a AT(TT/6)d 3 /l .33ndk a AT st (5-2) 

where d is the diameter of the spark kernel, AT s j is the temperature 
difference between the stoichiometric, adiabatic flame temperature of the 
kernel and the ambient spray and p a , c pa and k a are the density, specific 
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8.37E-02 

(1300*10 

2.58E-04 

( 1300'K) 


Table 5.1 


Fuel and Air Properties Used for Ignition 
Modeling 
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heat and thermal conductivity of air evaluated at 1300 K. This reduced to: 


t q - d 2 /6a 


(5-3) 


where a is the thermal diffusivity of air. 

The evaporation time was defined as the ratio of the mass of fuel 
contained within the spark kernel to the rate of fuel evaporation and was 
given by: 


t e * ( l-Q)p f c pa D 2 /8k a 01n( 1+B st ) (5-4) 

where 0 is the fraction of fuel In vapor form, D is the droplet diameter, 0 is 
the equivalence ratio and B s t is the stoichiometric combustion mass transfer 
number. The (1-0) term was included to account for fuel vapor in the 
ignition zone prior to the spark. This fuel vapor was assumed to have 
reduced the evaporation time by the time it wouid have taken to produce that 
vapor had it been present initially in the form of fuel droplets the same size 
as those present. Values of B S { were taken from Kanury ( 1970). 

The chemical reaction time was defined in terms of the laminar flame 
speed (S|_) and was given by: 


t c =12.5a/S L 2 (5-5) 

Using relations 5-3 to 5-5 in equation 5-1, the critical quenching distance 
was determined, and was given by: 


dg = [( l-O)p f O 2 /p a 01n( l+B st ) + ( 10a/S L ) 2 ] 0 - 5 (5-6) 


The minimum Ignition energy was defined as the energy required to raise 
the critical diameter kernel to the stoichiometric adiabatic flame 
temperature (T st ), and E mjn was given by: 

^min “ CpaPa^st^/SJdq^ (5-7) 


5.1.2 Minimum Ignition Energy Predictions 

Initial attempts to model n-heptane spray ignition were unsuccessful 
because the chemical times calculated with equation 5-5 were an order of 
magnitude higher than the evaporation times. This resulted in predicted 
ignition energies an order of magnitude higher than the experimental values. 
The chemical times were expected to be less than the evaporation times, 
since evaporation has been shown to be the rate controlling factor in ignition 
(Ballal and lefebvre, 1978; Peters and Mellor, 1981). For this reason, the 
chemical time (t c ) was neglected for further modeling attempts. (Values of 
the laminar flame speed (Warnatz, 1984) used in Equation 5-5 varied from 30 
to 40 cm/s, with maximum values occuring at 0 ® 1.1.) 

The spark ignition model, neglecting chemical times, was used to predict 
minimum ignition energies for the experimental spray cases reported 
previously. Figure 5.1 compares the predicted values of E m j n with the 
measured values for n-heptane and methanol sprays. Several observations 
can be made about the model from these results. The predicted effect of 
equivalence ratio matched the experimental results well for 0 > 0.8, as seen 
by the similarity of the slopes in this region. For equivalence ratios 
decreasing below 0.8, however, the model failed to predict the sharper 
increase in E m j n observed experimentally. This was expected since the 
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Figure 5. 1 Predictions of Minimum Ignition Energy using 
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model didn't Include chemical times, which would become larger as the lean 
Ignition limit was approached. For equivalence ratios above 0.8, the 
quantitative agreement between the predicted and experimental ignition 
energies was very good for the largest drop sizes. As droplet size 
decreased, however, the model increasingly underpredicted E m j n . This was 
partially explained by the exclusion of chemical times. As the droplet size 
decreased, the evaporation times decreased also and the chemical times 
became relatively more important. Furthermore, equation 5-6, without the 
chemical time, predicts a monotonic decrease in quenching distance and 
ignition energy, approaching zero for zero drop size (vapor). Therefore, 
the model also breaks down as droplet size decreases because it cannot 
account for the transition region fuel distribution effects which result in the 
optimum droplet size for ignition observed experimentally. 

In summary, the spark ignition model developed by.Ballal and Lefebvre 
was not able to account for spray ignition behavior as the equivalence ratio 
approached the lean limit or as the droplet diameter dropped below 40 urn. 
However, the general agreement between the model predictions and the 
experimental data is quite good considering the global nature of the model 
compared to the complex nature of the spark ignition process. 

5.1.3 Effect of Extent of Prevaporization 

The effect of fuel prevaporization on ignition was also studied using the 
ignition model of Ballal and Lefebvre. This was done to address a question 
raised in Chapter 4 concerning the fact that ignition energies were measured 
for n-heptane and methanol sprays with different values of Q. It was not 
clear how much of the difference in ignition energy between the n-heptane and 
methanol sprays could be attributed to their difference in 0, and how much 
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was due to the difference in their volatilities. To examine this, E m ^ n was 
predicted for sprays of the two fuels with the same value of Q = 0.39. This 
was done for n-heptane sprays with D = 41 urn and methanol sprays with D = 
44 ym. 

Figure 5.2 shows the effect of varying Q on the predicted ignition 
energies for these sprays. The experimental values of Q (n-heptane = 0.45; 
methanol = 0.32) are included as a reference. As seen, using a mean value 
of 0 (0.39) decreased the difference in E mjn between the n-heptane and 
methanol sprays by about 25%. This indicates that the majority of the 
difference in ignition energies was due to the difference in volatility of the 
fuels. 

5.2 Characteristic Time Model for Ignition (Peters and Mellor) 

The second ignition model used to evaluate the experimental data was the 
characteristic time model (CTM) for ignition of quiescent fuel sprays, which 
was developed by Peters and Mellor (1980) based on an earlier model for 
flame stabilization. This model associated specific times with the physical 
processes occuring during spark ignition and stated that the ignition limit W3S 
reached when the mixing rate of the spark kernel with its environment 
equalled the fuel evaporation rate. 

5.2.1 Model Formulation 

The basic assumptions of the model were similar to those presented for 
the previous model of Ballal and Lefebvre. A spherical spark kernel was 
formed by the ignition spark and the temperature in this kernel was the 
stoichiometric, adiabatic flame temperature of the fuel. The air properties 
in the kernel were evaluated at a mean temperature of 1300 K. In addition, 
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ignition was assumed, to be controlled by evaporation. That is, the time 
required for evaporation of the fuel in the spark kernel was very much larger 
than the time required for chemical reaction. Therefore, chemical times 
were neglected. The basic premise of the model was given by: 

t s i~t eb /0 (5-8) 

where t s j is the time required for the heat to be removed from the kernel by 
conduction, t eb is the time required for evaporation of a fuel droplet in the 
kernel and 0 is the equivalence ratio at the spark gap. Since the total amount 
of fuel being vaporized contributed to the ignition process, t eb was divided by 
the equivalence ratio, which was proportional to the total number of droplets 
present. 

Note that the CTM differed from the previously described model of Ballal 
and Lefebvre in that it only assumed proportionality between the relevant 
characteristic times, not absolute equality. The model was applied by 
evaluating t s j and t eb from the measured ignition energies and experimental 
conditions to determine the proportionality in equation 5-8. E m j n was then 
predicted based on this proportionality. The fuel and air properties used in 
these calculations were identical to those used in the previous model, and 
are listed in Table 5. 1 . 

For a quiescent mixture, the relevant mode of heat transfer from the 
spark kernel was conduction. Therefore, the time required for heat loss 
from the kernel was given by: 


*s! = A / a a 


(5-9) 


where A is the surface area of the kernel and a a is the thermal diffusivity. 
For a spherical spark kernel, this became: 

t S ] = TTd 2 p a c pa /k a (5 _ 10) 

where d is the diameter of the spark kernel and p a , c pa and k a are the 
density, specific heat and themal conductivity of air evaluated at 1300 K. To 
evaluate d, the minimum ignition energy was defined as the energy required 
to heat a spherical volume of air with diameter d to the stoichiometric, 
adiabatic flame temperature: 


dq = ^min/ TT/ ^Pa c pa AT st^ 1/ ^ (5- 1 1 ) 

where d has been replaced by dq, the quenching distance, and AT st is the 
temperature difference between the stoichiometric, adiabatic flame 
temperature of the kernel and the ambient spray. 

The time required for evaporation of a fuel droplet in the spark kernel 
was derived from the "d^ law" of Godsave, given by: 

<5-12) 

where D 0 is the droplet diameter and f) is the evaporation coefficient. This 
expanded to: 


*eb = ^o^Pf c pa/ 8l< a^ n ^ ,+ ®st^ ^ 13) 


where B s ^ is the stoichiometric mass transfer number, taken from Kanury 
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(1970). 

Using the experimental conditions and measured ignition energies, t s) and 
t eb were ca,cu,ate<1 with equations 5-10, 5-11 and 5-13. A simple linear 
regression analysis was then performed with these values to determine the 
proportionality in equation 5-8. Results of this are given in Figure 5.3, which 
shows t s j plotted as a function of t eb /0 for the experimental data of this 
study. The least-squares fit was: 


t s i = 15.8t eb /0 -0.7 (5-14) 

The coefficient of correlation (r) for this fit was 0.89 and the y-intercept 
was well within the standard error of estimate ( 13.6) of zero. The scatter in 
the correlation was not excessive considering only 26 data points were 
modeled. 

This correlation was very close to the one obtained by Peters and Mellor 
( 1980) using the CTM with experimental data of Ballal and Lefebvre ( 1978): 

t s j = 14.2t eb / 0 + 2.29 (5-15) 

Equation 5-15 was the result of Ignition energy measurements made for 
iso-octane, diesel oil and heavy fuel oil sprays with SMD ranging from 20 to 
170 ym and equivalence ratios ranging from 0.43 to 1.0. The agreement 
between the two correlations is very good considering t S ) and t eb only ranged 
from 0-120 ms and 0-7 ms respectively for the present study, while the 
corresponding ranges for the data in Peters and Mellor's correlation were 
0-1000 ms and 0-70 ms. 
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Figure 5.3 Characteristic Time Correlation for N-Heptane and 
Methanol Sprays (using Initial Droplet Diameters) 



5.2.2 Minimum Ignition Energy Predictions 

Using the correlation obtained for the present data (5-14) and neglecting 
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the small y-intercept, E m j n was calculated with equations 5-10, 5-11 and 
5-13 as: 

^min = °'261c pa p a ®'^AT s jO o ^(pf/01n( 1+B) ) ^ (5-16) 

Comparisons of the minimum Ignition energies predicted by 5-16 and the 
experimentally measured values are shown in Figure 5.4 for n-heptane and 
methanol sprays. As seen, the CTM qualitatively predicted the effects of 
droplet size and equivalence ratio on ignition quite accurately. Quantitatively 
however, the predicted ignition energies were about 50% higher than the 
measured values for n-heptane sprays of all sizes. For methanol sprays, 
the model overpredicted E m j n for the 57 ym sprays and underpredicted E mjn 
for the 30 ym sprays by about 30%, while the 44 ym predictions were very 
close to the experimental values. As in Ballal and Lefebvre’s model, the 
agreement between the measured and predicted values of E m j n deteriorated 
with decreasing equivalence ratio, suggesting the need for chemical kinetic 
times in modeling lean sprays. 

5.2.3 Effect of Initial Versus Soark Gap Drooiet Diameter 

The CTM was developed to model ignition energies based on initial spray 
conditions, since detailed information about spray parameters at the ignition 
location is not always available. (It was not clear whether CTM the 
correlation obtained by Peters and Mellor for Ballal and Lefebvre’s quiescent 
data (Eq. 5-15) was obtained using initial or spark gap diameters.) In this 
study, however, measurements of droplet diameter were made directly at 
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Figure 5.4 CTM Minimum Ignition Energy Predictions for 
N-Heptane and Methanol Sprays (using Initial 
Droplet Diameters) 
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the spark gap. Therefore, the droplet evaporation term in the CTM was 
modified to use these measured values, instead of the initial droplet 
diameters used previously. This gave a more realistic estimate of the time 
required for evaporation in the spark kernel. With this modification the 
droplet evaporation time became: 

t eb = DZ Pr < W 8k a l "< l +6 st ) (5-17) 

where D is the droplet diameter measured at the spark gap. The correlation 
between t S ) and t eb was determined as before and found to be: 

t sl = 2O.6t eb /0 + 1.7 (5-18) 

Figure 5.5 shows a plot of t S ] as a function of t eb /0 for this case. The 
correlation coefficient (r) of this least squares fit was 0.9 and the 
y-intercept was within the standard error of estimate (12.9) of zero. As 
seen in Figure 5.5, the evaporation times were decreased by about a factor of 
one-third compared to the initial diameter case (Figure 5.3), which 
increased the slope of the correlation by about 30%. 

Using this correlation, the minimum ignition energy was given by: 

E mln .O.389c pa p,-O-5AT st D o 8 (p f /0ln(l+B))'-5 (5-19) 

Values of E m j n predicted by this equation are compared to the experimentally 
measured values in Figure 5.6 for n-heptane and methanol sprays. As seen, 
this modification had no qualitative effect on the predicted ignition energies. 
Using the spark gap droplet diameter improved the n-heptane predictions but 
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Figure 5.6 CTM Minimum Ignition Energy Predictions for 
N-Heptane and Methanol Sprays (using Spark 
Gap Droplet Diameters) 




had very little effect on the methanol predictions. This was expected since 
n-heptane's greater volatility resulted in larger differences between the 
initial and spark gap droplet diameters than for methanol sprays. 

5.2.4 Effect of Extent of Prevaporization 

The CTM was further modified to include the effects of prevaporized fuel 
in the ignition zone prior to sparking. This fuel vapor was assumed to have 
reduced the evaporation time by the time it would have taken to produce that 
vapor had it been present initially in the form of fuel droplets (Ballal and 
Lefebvre, 1981). The droplet evaporation time became: 

t« b * < l-fl)t> 2 PfCpa/8 k aln< l+B st ) (5-20) 

where Q is the fraction of fuel in vapor phase and D is the droplet diameter 
measured at the spark gap. The correlation between t s j and t eb was 
determined as before and found to be: 

t $1 = 29.8t eb /0 + 5.5 . (5-21) 

Figure 5.7 shows a plot of t S | as a function of t eb /0 for this case. The 
correlation coefficient (r) of this least squares fit was 0.9 and the 
y-intercept was within the standard error of estimate (12.6) of zero. The 
evaporation times were decreased by about a factor of two compared to the 
initial diameter case, which nearly doubled the slope of the correlation. 

Using this correlation, the minimum ignition energy was given by: 


E min = O.676c pa p a -o * 5 AT st D 3 ((l-Q)p f /0ln(l+B)) K5 (5-22) 
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Figure 5.7 Characteristic Time Correlation for N-Heptane 
and Methanol Sprays (using Spark Gap Droplet 
Diameters and Extent of Prevaporization) 




Values of E m j n predicted by this equation are compared to the 
experimentally measured values in Figure 5.8 for n-heptane and methanol 
sprays. As seen, the qualitative agreement was again quite good with the 
exception of equivalence ratios approaching the lean limit. In addition, the 
quantitative agreement between the predicted and measured ignition energies 
improved significantly for the n-heptane sprays, while little change was seen 
for methanol. Again, this was expected due to n-heptane's greater 
volatility, which resulted in smaller drop sizes and more prevaporized fuel 
than methanol. 

In general, the CTM was shown to be quite reliable in predicting the 
ignition energies of the sprays examined in this study. It appears that use of 
local (spark gap) spray parameters, such as droplet size and extent of 
prevaporization, results in more accurate predictions compared to use of 
initial parameters. Inclusion of chemical times, while not attempted in this 
study, might also improve model performance for lean and smaller droplet 
size sprays. 
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Figure 5.8 CTM Minimum Ignition Energy Predictions for 

N-Heptane and Methanol Sprays (using Spark Gap 
Droplet Diameters and Extent of Prevaporization) 




5.3 Ignition Frequency Model 


The experimental ignition energy results illustrate the probabilistic 
nature of spray ignition. This was seen in the variation in ignition frequency 
with spark energy. An ignition frequency model (IFM) was developed to 
account for these probabilistic effects by assuming that the randomness of 
spray ignition is caused entirely by variations in two of the main parameters 
affecting ignition, the spark energy and the equivalence ratio in the spark 
gap. Random normal distributions were simulated for these parameters, and 
these distributions were used in conjunction with the previously described 
CTM for ignition (modified to use the spark gap droplet diameter and the 
extent of prevaporization) to predict the ignition frequency for a given spark 
energy and spray. The model was formulated as described below. 

The probability distribution functions of the equivalence ratio (number 
density) in the spark gap and the spark energy were quantified 
experimentally in terms of their standard deviations. The standard deviation 
of the number density measurements was 50%, and was determined by 
statistically analyzing number density measurements taken for a time period 
of 10 ms (the average time a typical droplet or group of droplets spent in the 
spark gap). The standard deviation of the spark energy measurements was 
10%. This information was used by a Monte Carlo random normal distribution 
simulator to generate arrays of 50 probable values of equivalence ratio (0j) 
and spark energy (E S p j) for a given mean equivalence ratio and mean spark 
energy case. E m j n ,i was calculated for each value of 0j, using the CTM, and 
compared to E S pj. If E m j n j was less than or equal to E S pj, successful 
ignition occured. If E m j n j was greater than E sp j, no ignition occured. The 
ignition frequency for that spark energy was the number of successful 
ignitions divided by fifty. This was repeated over a full range of spark 


energies to generate an ignition frequency curve for a given mean equivalence 
ratio spray. A complete listing of the API program used to run this program 
is given in Appendix C. 

Typical results of the model are given in Figure 5.9, which shows ignition 
frequency curves generated by the IFM compared to the ignition frequencies 
obtained experimentally for representative n-heptane and methanol sprays. 
As seen, the IFM predicted quite accurately the observed variation in ignition 
frequency with spark energy for these cases. 

Figure 5.10 shows IFM predictions for both rich (0 * 1.42) and lean (0 = 
0.44) 41 urn n-heptane sprays. Several things can be observed from these 
results. For the very lean case, the IFM substantially underpredicted the 
50% ignition energy, or E m j n . This was because the ignition frequency 
predictions were based on minimum ignition energies calculated with the CTM. 
Where the CTM underpredicted E m j n , such as lean cases approaching the 
ignition limit, the IFM overpredicted the ignition frequency for a given spark 
energy. The ignition frequency model, then, is limited by the accuracy of the 
ignition energy model (CTM) used. For the rich case in Figure 5. Iff; the 50% 
ignition energy predicted by the IFM was within 5% of the experimentally 
observed value. The slope of the predicted curve was less steep than the 
slope of the experimental curve, however. This issue will be addressed in 
the following paragraph. The ignition frequency modeling results shown in 
Figures 5.9 and 5.10 were representitive of results obtained for both fuels 
over the entire range of droplet sizes and equivalence ratios tested. 

The slope of the ignition frequency curves predicted by the IFM were a 
function of the standard deviations of the spark energy and equivalence ratio 
used in the model. Smaller standard deviations would have resulted in 
steeper slopes, and vice-versa. This is illustrated in Figure 5.11, which 
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shows the predicted Ignition frequency curves for a 41 um n-heptane spray 
with equivalence ratio of 1.42 using standard deviations of 50%, 30% and 10% 
for the equivalence ratio. The standard deviation of the spark energy was 
held constant at 10%. From Figure 5-11 it appears that by using a standard 
deviation of about 30% for the equivalence ratio in the IFM, the slope of the 
corresponding experimental ignition frequency curve in Figure 5.10 would 
have been predicted quite accurately. However, no decrease in the standard 
deviation of the number density with increasing equivalence ratio was 
observed experimentally to warrant this. 

Finally, note in Figure 5.11 that all three curves cross the 50% ignition 
line at approximately the same spark energy. This is because the probability 
distribution functions used in the IFM for both the equivalence ratio and the 
spark energy are symmetric about the mean value. If the standard 
deviations of the distributions were zero (i.e. if ignition was deterministic), 
the ignition frequency curve would collapse to a straight vertical line which 
would also cross the 50% ignition line at the same energy. This illustrates 
that while spark ignition of sprays is probabilistic, choosing 50% ignition as 
the criterion for the minimum ignition energy accounts for these probabilistic 
effects. 
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CHAPTER 6 


CONCLUSIONS AND RECOMMENDATIONS 


6. 1 Conclusions 

A monodlsperse spray Ignition system was developed and used to study 
the spark Ignition behavior of fuel sprays over a wide range of operating 
conditions. A criterion based on the probabilistic nature of ignition was 
determined for the minimum ignition energy, which was measured over a 
30-57 ym droplet range and a wide range of equivalence ratios for n-heptane 
and methanol fuels. The sprays were characterized in terms of their number 
density, drop size, extent of prevaporization and velocity at the ignition 
spark gap. The spray droplet size was an important factor in spray ignition, 
ignition was enhanced in all cases by decreasing the droplet diameter, which 
increased the total surface area of fuel in the ignition zone. Increasing the 
spray equivalence ratio also enhanced ignition in all cases, so the optimum 
amount of fuel In the ignition zone was never reached. Fuel volatility was the 
critical fuel property in correlating spray ignition behavior for a given fuel 
droplet size. 

Minimum ignition energy measurements were made for prevaporized, 
premixed n-heptane and methanol to provide ignition data for the lower limit 
of droplet size. The optimum vapor phase equivalence ratio for ignition was 
also determined from these measurements for n-heptane. Activation energy 
was the critical fuel property for correlating prevaporized ignition behavior. 

Flame propagation was closely related to ignition. The importance of 
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transition region effects on flame propagation rates in sprays of this size 
range and smaller was exhibited by the existence of an optimum droplet size 
for ignition at lean equivalence ratios for both fuels. This optimum size 
could not be quantified due to limitations of the spray generator, but it was 
attributed to the maximum flame speeds achieved by sprays in the transition 
region. Extension of the lean ignition limits was observed for sprays of both 
fuels, and attributed to the enhancing effect of the fuel vapor supplied by 
evaporating droplets. 

Two models using a characteristic time approach were applied to 
correlate the spray ignition results. Both models predicted quite 
accurately, qualitatively and quantitatively, the observed ignition behavior. 
The models underpredicted ignition energies as the spray equivalence ratios 
approached the lean limit and as the droplet sizes decreased below 40 urn. 
This was because they did not account for the importance of kinetic times for 
lean mixtures and smaller droplets, nor did they account for transition 
region effects in smaller sprays. 

The random nature of the spray ignition process was explained with a 
probabilistic model which took into account the distributions of spray number 
densities and spark energies during ignition. This model predicted ignition 
frequency curves which were in very good agreement with the observed 
variation in ignition frequency with spark energy. The success of the model 
indicated that the probabilistic nature of spray ignition is due in a large part 
to the random positioning of droplets in a spray. 



6.2 Recommendations for Future Work 


While this study points out important factors for optimizing ignition and 
minimizing ignition energies for spray combustors, further experimental and 
theoretical advancements are needed before the knowledge acquired here can 
be applied to such systems. Suggested experimental studies and their 
objectives are listed below: 

1) Modification of the spray system should be done to generate droplets in 
the 10-30 urn range to quantify the optimum drop size for ignition and to 
determine the effects of equivalence ratio, prevaporization and fuel 
properties on this optimum. 

2) Ignition studies should be performed with the prevaporizing system to 
add vapor to the interdroplet spacing of the sprays. 0 should be varied for 
fixed values of D and 0 to quantify the importance of prevaporization to 
ignition and to evaluate ignition models which use Q. 

3) Ignition studies should be performed with sprays formed using the 
dispersion jet only (no dilution air) to reduce the levels of turbulence in the 
sprays. A larger diameter test section should also be used to eliminate the 
problem of droplet impingement on the walls. 

4) The prevaporizing apparatus should be modified to generate rich 
methanol mixtures in order to quantify the optimum equivalence ratio for 
prevaporized methanol mixtures. This could possibly be done by heating the 
entire flow system from the vaporizing chamber to the test section, or by 


increasing the fuel and air flow rates for the existing set-up. 

5) Modifications should be made with the spray system to generate bi-modal 
distributions of droplets to quantify their ignition behavior relative to 
monodisperse sprays. Such a study would begin to provide a link between 
monodisperse and polydisperse ignition studies. This information could also 
be used to assess the relative importance of size correlation parameters, 
such as SMO and Dzq, for ignition modelling. 

Improvements in modeling of the spark ignition process can be achieved 
in several ways: 

1 ) A more detailed chemical kinetic term is needed to account for the sharp 
increase in ignition energies observed experimentally as the lean ignition 
limit is approached. The use of a non-unity exponent for 0 in the t^ c term 
developed by Peters and Mellor ( 1982) could improve the CTM performance in 
this respect. 

2) Attention should be given to modelling the observed transition region 
effects, i.e. the flame propagation and fuel distribution effects resulting in 
an optimum droplet size for ignition. Both ignition models presently predict a 
monotonic decrease in ignition energy with decreasing droplet size, thus 
underpredicting E m j n for smaller droplets (D < 35 urn). 

3) Statistical methods for determining the actual distributions of number 
density and equivalence ratio for sprays over a range of conditions should be 
developed in order to verify the IFM assumptions in this respect. 
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APPENDIX A 


COMPREHENSIVE SPRAY CHARACTERIZATION AND IGNITION RESULTS 

For brevity, only representitive results of the spray characterization 
and ignition measurements were reported in the main body of text. 
Comprehensive results are given in this appendix. Figures A. 1 and A. 2 show 
radial profiles of the local equivalence ratio for all of the sprays tested in 
this study. Figure A. 3 shows Fraunhofer diffraction patterns from each of 
the fuel and droplet size sprays tested at an overall equivalence ratio of 0.7. 
Finally, Figures A. 4 and A. 5 report the experimental ignition energy 
measurements, showing the variation in ignition frequency with spark energy 
for each of the sprays tested. 
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Figure A. 1 Radial Profiles of Local Equivalence Ratio for 

N-Hept.ane Sprays of Varying Overall Equivalence 
Ratio (D 0 = 65, 50 and 40 ym) 
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Figure A. 2 Radial Profiles of Local Equivalence Ratio for 

Methanol Sprays of Varying Overall Equivalence 
Ratio (D 0 = 65, 50 and 35 jam) 


& 




ORIGINAL pA 

OF POOR QUALITy 122 




H — 53 



H — 41 


M — 57 



M — 44 


H — 33 


M — 30 




Figure A. 3 Fraunhofer Diffraction Patterns produced by 

N-Heptane (D = 53, 41 and 33 um) and Methanol 
(D = 57, 44 and 30 um) Sprays (0 O = 0.7) 
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Figure A. 5 Effect of Spark Energy and Equivalence Ratio 
on the Ignition Frequency of Methanol Sprays 
(D = 57, 44 and 30 um) 





APPENDIX B 


COMPUTER PROGRAM F( 


)ATA CQLI 


)N AND ANALYSIS 


PROGRAM LDA4 

Q********1Ht******************************************** 

C EXPONENT = 8; NCYCLES = 2" 4 
C THIS PROGRAM ACQUIRES AND ANALYZES DATA FROM A 
C TSI MDL 1 990 LDV SYSTEM ; THE LDV DATA READY PULSE 
C IS USED TO TRIGGER CH.O A/D VELOCITY DATA CONVERSION 
C THROUGH THE A/D EXTERNAL TRIGGER ; THE CLOCK COUNTS 
C THE TIME FOR NSAMP DATA, GIVING THE DROPLET RATE; 

C THE LDV SETTINGS ARE: CYCLES/BURST=2"4; EXP=8; 

C COMPf =7 ; FREQ= 1 00: 1 ; GAIN= 1 ; AMP LIMOFF; 

C FILTERS: HI=3MHZ; L0=30KHZ; 

0 ************#**»***#*********»*»*******#****-*****#*** 
INTEGER*2 I BELL 
DATA IBELL/"07/ 

DIMENSION IDAT( IOO),DAT( IOO),YOLT( 100),VEL( 100) 

DIMENSION DFREQ( 100),VELMN( 100),YELRMS( 100),RM$P( 100) 

Q***************************************************** 

C MACRO DATA TRANSFER ROUTI NE ( ADCLS.MAC) ARGUMENTS: 

C I DAT: ARRAY FOR VELOCI TY DATA 
C ICNT: VARIABLE FOR FREQUENCY DATA 

C NSAMP: NUMBER OF DATA TAKEN ( UP TO 1 000) 

Q*«**«******«****«*******«*******************«****«*** 

WRITE(5,96) 

96 FORMATOOENTER THRESHOLD VELOCITY ( M/S):‘) 

READ(5 ,*)VTH 

97 CONTINUE 
NSAMP=100 
WRITE(5,98) 

98 FORMAT( ’ ENTER NAVG ( I - 1 00 ):' ) 

READ(5,*)NMN 

99 CONTINUE 
VBUF=0. 

FBUF*0. 

RBUF=0. 

DO 103J-I.NMN 

CALL ADCLS( I DAT ,ICNT .NSAMP ) 


CALCULATE MEAN VELOCITY (M/S) 
NOTE: VEL=FREQ*FRINGE SPACING 
FREQ=YOLT*( 16*10"! 0/2"( n+ 1 0)) 
n=8 ; FREQ=V0LT*6 1 0352 
FRINGE SPACING=3.3 MICRONS 
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I 

I 


SUM=0. 

DO IOOM.NSAMP 

YOLT( I )=5.*FL0AT( IDAT( I ))/4095. 

VEL( I )=VOLT( 0*2.0142 
IF (VEL(O.OT.YTH) GO TO 100 
SUM=SUM+VEL(I) 

100 CONTINUE 
VELMN(J)=SUM/NSAMP 

£*■)«***»**************«*********«*******«****«****** 
C CALCULATE RMS VELOCI TY ( M/S) 

q»****************«************ ********* ************ 

SUMM=0. 

DO 101 1=1 ,NSAMP 

SUMM=$UMM+((VEL(I)-VELMN(J))**2) 

101 CONTINUE 

SUMMM=SUMM/FLQAT(NSAMP) 

VELRMS(J)=SQRT(SUMMM) 

RMSP(J)=VELRMS(J)/VELMN(J) 

Q*************************************************** 

C CALCULATE DROPLET RATE (KHZ): (CLOCK RATE=0. 1 KHZ) 

Q*************************************************** 

RATE=0. 1 

DFREQ(J)=FLOAT(NSAMP)/(FLOAT(ICNT)/RATE) 

Q*****************************«********************* 

C STORE VALUES FOR MEAN 

Q*************************************************** 

YBUF=YBUF +YELMN( J) 

FBUF=FBUF+DFREQ(J) 

RBUF=RBUF +DFREQ( J) 

103 CONTINUE 

Q*************************************************** 

C CALCULATE MEAN VALUES 

Q*************************************************** 

VELMNM=VBUF/FLOAT( NMN) 

DFREQM=FBUF/FLOAT( NMN) 

TINTM=RBUF/FLOAT(NMN) 

WRITE( 5 , 1 04)DFREQM .YELMNM .TINTM 

104 FORMATC MEAN YALUES (N,V ,T): ‘ .3F7.3) 

WRITE(5,1000)IBELL 

1000 FORMATC+'^l) 

GO TO 97 

STOP 

END 


SUBROUTINE ADCLS.MAC: 


AOCLS: 


LOOP: 


.TITLE ADCLS 
.GLOBL ADCLS 
ADCSR* 1 70400 
ADBUF* 170402 
CCSR= 1 70420 
CBUF* 170422 
CLRADCSR 
CLRADBUF 
CLRCCSR 
CLR CBUF 
MOV 2(RS),RI 
MOV 4(R5),R2 
MOV ®6(R5),R3 
BIS*55,CCSR 
BIS*20ADCSR 
TSTB ADCSR 
BGE LOOP 

MOYADBUF.(RI)* 
SOB R3, LOOP 
BIS*1000,CCSR 
MOV CBUF ,(R2) 
CLRADCSR 
CLRADBUF 
CLRCCSR 
CLR CBUF 
RTS PC 
.END 


ADCLS.MAC MAKES 1 000 A/D 
.•CONVERSIONS ON CH.O 
,-OF VELOCITY DATA FROM 
;THE ANALOG OUTPUT OF 
;THE LDV TIMER. IT ALSO 
.RECORDS THE * OF COUNTER 
;TICS DURING ACQUISITION 
;T0 GIVE DROPLET FREQUENCY 


IVOLT ARRAY LOCATION 
ICNT VARIABLE LOCATION 
SET NSAMP 

SET RATE=0. 1KHZ, M0DE=2, GO 
SET A/DCH=0, EXT. TRIG. ENABLE 
CHECK IF CONV. DONE 
LOOP TILL DONE (NEGATIVE) 

PUT IVOLT(I) INTO ARRAY 
COUNT* OF CONVERSIONS 
SIMULATE ST2 INPUT 
READ COUNTER* INTO ICNT 



APPENDIX C 


APL IGNITION FREQUENCY MODEL PROGRAM 


<7 1 GN I T I OH ;N ;PH I BfiR ;PH I STD jPH I ;E11 1 N;DROP ; OMEGA ;CPA ;RHOfl ;LM8 ;DTST ;PROB ;COMB ;M 
IIJ ’ INPUT THE SPARK GRP DROPLET DIAMETER = * 

(2 1 DROP - 0 
131 ' INPUT OMEGA- ' 

14] OMEGA - □ 

153 ' INPUT THE MEAN SPARK GAP EQU I UALENCE RATIO* • 

16] PHIBAR - □ 

17] 'INPUT THE STANOARO DEU I AT I ON OF THE SPARK EQU I UALENCE RATIO (PERCENT) : ' 
181 PH I STD - □ 

19) N - 50 
MO] M - NX28 

MU 'INPUT THE STANDARD DEU I AT I ON OF THE SPARK ENERGV (PERCENT)*' 

I 121 ESTD - 0 
I 13] PHI - SIMNORM M 

M4] PHI - PHIBAR + (((PHISTOx.ODxPHIBAR)xPHI ) 

1151 PHIT - PHIx(PH1>0) 

M6J PHI - PHIT + ((PHIT=0)x. 01) 

[17] PHI - (28 50)PFHI 
MSI ESPARK-SPARKV ESTD 
I 19] CPA - 1200 
1201 RHOA * .271 
(211 DTST - 1983 
[22] RHOF - 688 
1231 LNB - 1.92 

124] EMIN - .676xCPAxDTSTx(RHOA*-0.5) 

I25J EMIN - EH I Nx ( ( RHOF+LNB ) * 1 . 5 ) x ( 0ROP*3 > x ( < ( 1 -OMEGA ) +PH I ) * 1 . 5 ) 

[26] EMIN - EM I NX IE" 15 

[27] PROB - (+/(EMIN<ESPARKmN 
128] PROB - <28 DpPROB 

1291 " 

[301 

1311 'THE SPARK ENERGIES AND IGNITION PROBABILITIES ARE*' 

1321 ESPMN, PROB 

133] COMB - ESPHN, PROB 

134] NUMPUT COMB 

c2 - SFfiRKV ESTD;ESPEXP;ONE;H;ESFARK 
l U ONE - <1 50) PONE*- 1 
[2 J ESPEXP - ESPMN +.x ONE 

13] M «- 50x28 

141 ESPARK - SIMNORM M 

[51 ESPARK «- (28 50)PESFARK 

[61 2 *- ESPEXP + (<<ESTDx.01)xESPEXP)xESPRRK) 

<?Z*-S I MNOR11 N;A;B;THETA;R 
[ 1 J N-N+2 
[21 fT-UN I FORM N 
131 B*-UN I FORM N 

14] THETA- < 02 )xA 

15] R-C2x(«B))*0.5 

161 2-(Rx( loTHETA) ) , <Rx< 20 THETA) ) 
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